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ABSTRACT 
In  this  report,  the  basic  heat  transfer  and  isothermal  flow  friction 
characteristics  for  five  perforated  nickel  matrices,  one  perforated  nickel 
fin  with  solid  splitter  matrix,  and  one  solid  nickel  matrix,  all  of 
"parallel  plate"  configuration  are  presented.   In  addition,  similar  in- 
formation is  provided  for  a  20°  skew  matrix  fabricated  from  perforated 
nickel.   Finally,  performance  of  two  offset  rectangular  fin  matrices 
constructed  from  stainless  steel,  is  evaluated.  Where  meaningful,  a 
comparison  of  the  performance  of  these  surfaces  is  made.  All  matrices 
are  investigated  by  the  transient  method  of  testing  utilizing  the  maximum 
slope  analysis  of  Locke  and  include  the  effects  of  longitudinal  thermal 
conduction  as  reported  by  Howard. 
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NOMENCLATURE 


English  Letter  Symbols 

2 
A        Matrix  total  heat  transfer  area,  ft 

2 

A        Matrix  minimum  free  flow  area,  ft 

2 
Afr       Matrix  total  frontal  area,  ft 

A8        Solid  matrix  cross  sectional  area  available  for  thermal 
conduction,  ft 

2 
A.        Conduction  area  corrected  for  effect  of  perforations,  ft 

A*  Plane  surface  area,  ft^ 

a  Plate  thickness,  ft 

a  Short  side  of  a  rectangular  flow  passage,  ft 

b  Plate  spacing,  ft  or  in. 

b  Long  side  of  a  rectangular  flow  passage,  ft 

C  Flow  stream  capacity  rate,  (mc  ),  Btu  (hr  °F) 

(T  Matrix  capacity,  W  c  ,  Btu/  °F 

o  S  S 

c         Specific  heat  at  constant  pressure,  Btu/ (lb  °F) 
P 

c         Matrix  material  specific  heat,  Btu/(lb  °F) 
D„        Hydraulic  diameter,  4r^ 

d         Inside  pipe  diameter,  in. 

2 
E         Flow  friction  power  per  unit  area,  HP/ft 

f         Mean  friction  factor,  dimensionless.   This  is  the  "small" 
or  "Fanning"  friction  factor.   (Ratio  of  wall  shear  stress 
to  the  fluid  dynamic  head.) 

G         Exchanger  flow  stream  mass  velocity,  (m/A  ),  lb/(hr  ft*) 

c 

g         Proportionality  factor  in  Newton's  Second  Law,  g  =  32.2 
c        (lb/#)  (ft/sec2)  C 

h         Unit  conductance  for  thermal  convection  heat  transfer, 
Btu/(hr  ft*  °F),  or  heat  transfer  power  per  unit  area 
per  degree  temperature  difference,  Btu/(hr  ft^  °F) 


ix 


English  Letter  Symbols  (continued) 

2/3 
j        Colburn  factor  =  N   NpB    ,  heat  transfer  characteristic, 

dimensionless 

j'        Colburn  j  factor  based  on  plane  surface  area,  A*,  dimension- 
less 

K  , K       Contraction  loss  coefficient  for  flow  at  heat  exchanger 
entrance  or  exit  respectively,  dimensionless 

k         Unit  thermal  conductivity,  Btu/(hr  ft  °F/ft) 

kg       Matrix  thermal  conductivity,  Btu/(hr  ft  °F/ft) 

L        Total  matrix  flow  length,  ft 

m        Mass  flow  rate,  lb/hr 

2 
P        Pressure,  #/£t 

p  Porosity  for  matrix  surfaces,  dimensionless 

q  Heat  transfer  rate,  Btu/hr 

R  Universal  gas  constant,  ft-#/lb-  ° R,  (53.35  for  air) 

r.  Hydraulic  radius,  (A  L/A),  ft,  (4r^  =  hydraulic  diameter) 

s  Solidity  of  a  perforated  material,  1-3^  dimensionless 

T  Absolute  temperature,  degrees  Rankine,  °R 

t  Temperature,  degrees  Fahrenheit,  °F 

V„        Matrix  volume,  ftJ 

m  ' 

V        Material  volume  corrected  for  effects  of  perforations 
s  r 

W        Mass  of  matrix,  lb 
s 

x         Distance  along  the  flow  passage  in  direction  of  flow,  ft 

Greek  Letter  Symbols 

°^        Aspect  ratio  of  a  rectangular  flow  passage,  b/a,  dimension- 
less 


0 


Compactness;  ratio  of  total  heat  transfer  area  to  the  volume 
between  the  plates,  ftVft 


A       Denotes  difference 
°        Fin  thickness,  ft 

x 


Greek  Letter  Symbols  (continued) 

J  Conduction  area  reduction  ratio  due  to  perforations  in 

perforated  material;  cross  sectional  area  of  perforations/ 
cross  sectional  area,  d imensionless 


0 


Time 


A>       Longitudinal  conduction  parameter,  k  A  ,  d imensionless 

s  s 

mLc 
P 

1       Time  parameter 

JUL      Fluid  viscosity  coefficient,  lb/hr-ft 

~y  Area  reduction  ratio  due  to  perforations  in  perforated 

material,  dimensionless 

p      Density,  lb/ft3 

O"      Ratio  of  freeflow  area  to  frontal  area,  Ac/Afr,  dimension- 
less 

P  Denotes  "function  of" 

A      Compactness  for  perforated  material  including  tbe  effect 

of  area  reduction,  A*s/V  .  dimensionless 

*  m' 

/_.  Summation 
Subscripts 

atm  Local  atmosphere 

f  Fluid 

i  Initial,  individual 

k  Equivalent 

m  Mean  or  matrix,  as  appropriate 

0  At  orifice 

s  (Matrix  material)  solid 

STD  Standard  temperature  and  pressure 

1  Upstream 

2  Downstream 
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Force  and  Mass  Units 

lb       Denotes  pounds  mass  to  differentiate  from 

#        Denoting  pounds  force 
Dimensionless  Groupings 

N_        Reynolds  number,  (4r  G/a  ),  a  flow  modulus 

Np,  v     Reynolds  number  for  pipe,  (dG/yu  ) 

N_fc       Stanton  number,  (h/Gc  ),  a  heat  transfer  modulus 
St  p 

NNu       Nusselt  number,  (h4r  /k),  a  heat  transfer  modulus 

Np        Prandtl  number,  (uc./k),  a  fluid  properties  modulus 

N-        Number  of  heat  transfer  units,  (hA/mc  ) 

j         Generalized  heat  transfer  grouping  -  Colburn  "j"  factor, 

(N.  N  2/3).   This  factor  versus  NR  defines  heat  transfer 
characteristics  of  the  surface. 

f        Mean  friction  factor.   This  is  the  "small"  or  "Fanning" 
friction  factor  (Ratio  of  wall  shear  stress  to  the  fluid 
dynamic  head).   This  factor  versus  N  defines  the  fric- 
tion characteristics  of  the  surface. 


A.        Conduction  parameter,  k  A  /mLcp   for  solid  material; 

P 


r 


k_A. /mLc_  for  perforated  material 


\         Equivalent  conduction  parameter  (corrected  for  equivalent 
length  in  perforated  mater ial).x  L 


Time  parameter,    hA  q 

W  c 
s  s 
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1.    INTRODUCTION 

While  it  is  appreciated  that  the  designer  of  a  compact  heat  exchanger 
is  confronted  with  many  qualitative  evaluations,  nevertheless  he  is  general- 
ly working  specifically  within  limitations  relating  to  flow  frontal  area, 
weight,  pressure  drop  and  volume.   It  is  for  this  reason  that  different 
heat  transfer  surfaces  may  be  effectively  compared  on  graphs  using  co- 
ordinates of  f  and  j  (or  f/j  ratio)  versus  ND. 

As  in  many  other  disciplines,  model  theory  is  the  most  satisfactory 
approach  from  a  cost  effectiveness  standpoint.   The  use  of  such  methods 
as  steady  state  "steam  to  air"  tests  or  "transient"  testing  techniques 
provides  the  basic  flow  friction  and  heat  transfer  characteristics  without 
necessitating  the  construction  of  a  full  scale  heat  exchanger. 

The  data  presented  herein  was  obtained  by  the  use  of  a  transient  test 
technique  or  "single  blow".  C7l  This  technique  is  one  of  heating  a  sample 
core  to  some  predetermined  value  and  then  after  thermal  equilibrium  is 
achieved  in  the  specimen,  instantaneously  removing  the  upstream  heat  source 
and  monitoring  the  downstream  temperature  response  of  the  fluid  at  the  core 
exit. 

The  objectives  of  this  thesis  specifically  were:   to  modify  an  experi- 
mental test  rig  to  operate  with  a  larger  capacity  blower  and  to  test  sever- 
al offset  rectangular  fin,  of  stainless  steel,  and  "modified  parallel"  plate, 
of  perforated  nickel,  sample  cores. 

In  the  testing  of  perforated  nickel  specimens,   the  effect  of  varying 
shape  and  size  of  perforations  was  to  be  investigated. 


2.    SUMMARY  OF  THEORY 

The  "starting  stage  of  operation11,  "single  blow",  or  "transient  test 
technique",  as  it  has  been  called  by  various  authors,  had  its  inception 
with  the  analytical  treatment  by  Anzelius  in  1926.   Nusselt  in  1927,  Hausen 
in  1927  and  1929,  and  finally  Schumann  in  1929  treated  the  problem  of  ana- 
lytically determining  the  fluid  and  solid  temperature  as  a  function  of 
position  and  time  for  a  porous  solid,  initially  at  a  uniform  temperature, 
when  subjected  to  a  change  in  the  entering  fluid  temperature.   Some  repeti- 
tion occurred  as  the  authors  were  unaware  of  the  work  performed  by  their 
contemporaries  at  the  time  C9]  .   Anzelius  and  Schumann  made  similar  as- 
sumptions; however,  the  results  are  different.   Anzelius  simplified  the 
problem;  whereas,  Schumann  proceeded  to  obtain  an  explicit  mathematical 
solution  for  the  temperature  distribution.   This  solution  was  for  a  liquid 
flowing  through  a  porous  solid  and  consisted  of  two  infinite  series  which 
were  related  to  Bessel  functions.   Schumann  states  C173  that  if  the  fluid 
is  a  gas,  the  solution  to  the  problem  is  more  formidable.   He  proceeds  to 
relate  how  dimensional  analysis  may  be  used  so  that  his  investigation  re- 
mains usable  for  compressible  fluids. 

The  porous  solid  is  considered  to  be  homogeneous  and  at  a  uniform 
temperature  throughout  with  the  fluid  flowing  at  this  same  temperature. 
With  a  sudden  change  in  the  entering  fluid  temperature,  the  temperature 
of  the  fluid  and  the  solid  as  a  function  of  both  position  and  time  re- 
mains to  be  evaluated.   The  theoretical  generalized  heating  curves  are 
developed  from  an  energy  balance  on  the  porous  solid  as  presented  in  [19j 
and  restated  in  Cll  .   Inasmuch  as  experimental  procedure  used  in  this 
report  is  dependent  upon  the  assumptions  and  boundary  conditions  of  this 
energy  balance,  it  is  briefly  summarized  below: 
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ro 


CP"^ 


hb(trts)o(x 


* 


-ksAs|ks    c£> 


dx 


K\50LID 


ncp(v^<»«) 


m 


S 


^>  Wfc*ft)* 


Assumptions  made  in  the  analysis  are: 

a.  Properties  of  the  fluid  are  temperature  independent 

b.  Fluid  flow  is  steady 

c.  Porous  solid  is  homogeneous 

d.  Thermal  conductivity  of  both  fluid  and  solid  is  infinite 
perpendicular  to  the  fluid  flow  direction 

e.  Thermal  conductivity  of  fluid  is  zero  in  the  flow  direction 
Boundary  conditions  are  Cl33  : 

(a)  The  matrix  is  initially  at  a  uniform  temperature 

(b)  At  time  equal  zero,  the  temperature  of  the  entering  fluid 
changes  instantaneously  to  a  different,  constant  value, 
i.e.,  a  step  change  in  fluid  temperature 

(c)  The  matrix  boundaries  are  adiabatic 
Heat  absorbed  by  solid  =  Ps^scs  ^— 5.  dx 

Heat  transferred  to  the  solid  by  convection  =  nb  vtf"*-S' 

Heat  transferred  from  the  fluid  by  convection  =  mCa^P  ci* 

dx 

Heat  transferredfrom  solid  by  conduction  ■  -ksAvs  d±s  dx 

dx2- 

Thus  an  energy  balance  for  the  fluid  and  solid  yields: 


(1)  ^AsCs^dx  =:  ksAs^s  dx  +  hb(tf--ts)^* 

(2)  mcp|tfdx  +  hb(-tf--ts)dx  =  0 


Rearranging  the  above  equations  and  substituting  in  the  following  para- 
meters [3J  : 


T  =  generalized  time  variable  =  ■****• 

Wsc5 


e-  ^f 

m  L 


where  h  =  unit  conductance  for  convective  heat  transfer, 

Btu/hr  ft2  °F 

.2 


A  =  matrix  heat  transfer  area,  ft' 

Wgcg  =  matrix  capacity,  Btu/  °F 

"  =  time,  hrs. 

W  =  mass  of  fluid  contained  within  the  matrix,  lb. 

m  =  mass  flow  rate  of  fluid,  lb/hr. 

x  =  distance  from  entrance  of  matrix  along  flow  passage 
in  the  direction  of  flow,  ft. 

L  =  matrix  flow  length,  ft. 

Regrouping  the  above  equation  for  T"  and  introducing  the  specific  heat  of 

the  fluid,  cp,  we  get:  ^  _  ^^    _  h^WfCp  X 

WSCS     m  WSCSC?L 

Since  the  thermal  capacity  of  the  fluid  contained  at  any  time  within  the 
matrix  (Wfc  )  is  much  less  than  the  thermal  capacity  of  the  matrix  (W  c  ), 
the  second  term  can  be  neglected;  accordingly: 


Wscs 


9 


(3) 


z  =  generalized  position  variable  =  b O  — 

™Cp  L 

Inasmuch  as  the  number  of  transfer  units  (NT  )  =  J?   -      ,  this  simply 
becomes 

■■  NTU-*-  (*) 


in  c  L 
P 

2 
where  k  ■  matrix  thermal  conductivity,  Btu/hr  ft  °F  /  ft. 

A  =  solid  matrix  cross  sectional  area  available  for  thermal  conduc- 
s   tion,  ft2. 

m,  c  ,  and  L  are  as  given  above. 
P 

Equations  (1)  and  (2),  together  with  (3),  (4),  and  (5)  yield: 


(6) 


--£  =  "ts""^  (7) 


where  t  =  temperature  of  the  solid 
tf  -  temperature  of  the  fluid 
For  the  case  of  thermal  conduction  in  the  solid  parallel  to  the  direction 
of  flow  equal  to  zero,  equations  (6)  and  (7)  simplify: 


&-v±. 


i^  =  -ts--tf  w 


It  is  with  these  two  differential  equations,  that  Schumann  utilizing 

Bessel  functions  determined  the  time  -  temperature  history  for  the  solid 

after  a  step  change  in  fluid  temperature.   For  a  more  thorough  treatise, 

the  reader  is  referred  to  references  [13}  or  [17]  . 

The  task  of  determining  the  N   value  for  a  test  point  required  the 

Tu 

superposition  of  the  experimentally  obtained  time  -  temperature  heating 


or  cooling  curve  on  the  family  of  theoretical  curves.  Although  this  would 
uniquely  determine  the  N  value,  it  can  be  appreciated  as  being  a  tedious 
effort. 

Fortunately,  in  1950,  Locke  [13]  developed  a  "maximum  slope"  technique 
which  bypasses  this  task.   He  differentiated  the  theoretical  curves  and 
plotted  the  slope  of  the  curve  versus  T  /N   with  N   as  a  parameter.  The 


expression  for  the  slope  is 


d[VNTa]   l^r' 


-T.^vfFMrvTje 


-(NwV) 


L. 
where  the  downstream  fluid  temperature  is  evaluated  at  x  «■  L  so  that  z  = 

NL  and  t.  -  t.-. 
Tu      f    f2 

In  the  efforts  of  both  Locke  and  Schumann,  the  thermal  conductivity 
of  the  solid  parallel  to  the  flow  was  taken  as  zero.   Although  the  effect 
of  longitudinal  conduction  is  not  pronounced  at  high  Reynolds  number  range, 
it  is  considerable  at  low  laminar  Reynolds  numbers.   Several  authors  have 
considered  the  effects  of  longitudinal  conduction,  introducing  the  conduc- 
tion parameter,  A-  •   When  A.t  O  t   the  maximum  slope  is  not  a  unique 
function  of  N   but  depends  on  A  as  well.   Howard  [7]  has  tabulated  maxi- 
mum slopes  as  functions  of  N   and  A.  ,  as  well  as  presenting  graphs, 
which  he  evaluated  by  a  finite  difference  technique  with  the  aid  of  a 
digital  computer   The  evaluation  of  N   in  this  report  is  based  on  these 
tabulated  values  in  so  far  as  practicable. 

Satisfying  the  aforementioned  idealizations,  initial  and  boundary 
conditions,  leaves  the  determination  of  the  conduction  parameter,  A-   , 
and  the  experimental  evaluation  of  the  maximum  slope  of  the  heating  or 
cooling  curve  as  the  only  two  remaining  items  necessary  to  determine  the 
heat  transfer  characteristics  of  the  matrix. 


3.    EXPERIMENTAL  TECHNIQUES 

The  results  of  Howard's  finite  difference  calculations  in  reference 

[73  were  used  in  this  report  for  evaluating  N   and  ultimately  the 

Tu 

Colburn  j  factor.   For  this  reason,  the  idealizations  imposed  by  Howard 
were  restrictions  placed  on  the  experimental  set  up  used.   Howard  £7] 
considered  transient  heating  of  a  porous  solid  when  subjected  to  a  step 
change  in  fluid  temperature  including  the  material  longitudinal  conduc- 
tion effects.   The  idealizations  are  restated  below: 

"(1)  The  flow  entering  the  matrix  is  steady  and  uniform  in 

velocity  and  temperature,  and  remains  steady  and  uniform 
at  any  cross-section  as  it  flows  through  the  matrix.   The 
thermal  conductivity  of  the  matrix  is  infinite  in  the 
direction  normal  to  the  flow  and  finite  in  the  direction 
parallel  to  the  flow.   Thus  the  problem  is  one-dimensional 
in  space. 

(2)  The  thermal  capacity  of  the  fluid  contained  at  anytime 
within  the  matrix  is  small  compared  with  the  thermal 
capacity  of  the  matrix.   That  is,  the  fluid  is  normally 
restricted  to  a  gas  and  the  results  to  be  presented 

.will  not  be  expected  to  apply  for  liquid  fluids.   This 
also  means  that  for  the  fluid  there  will  be  no  time  de- 
pendent terms  in  the  equations. 

(3)  The  thermal  properties  of  the  fluid  and  matrix  are  con- 
stant and  uniform. 

(4)  The  convective  heat  transfer  coefficient  is  some  suitable 
average  and  remains  uniform  and  constant. 

(5)  At  time  zero  the  change  in  fluid  temperature  will  be  a 
step  function,  with  the  matrix  and  its  entrained  fluid 
initially  at  some  uniform  and  equal  temperature." 

Idealization  (1)  was  achieved,  as  far  as  steadiness  and  uniformity  of 

temperature  and  velocity  were  concerned,  by  responsible  experimentation. 

In  an  effort  to  improve  on  the  previously  used  test  apparatus  of  Bannon 

[l^  ,  a  new  inlet  cone  was  designed  and  the  test  section  relocated.   The 

inlet  cone  was  fashioned  according  to  the  specifications  of  Smith  and 


Wang   L18j  for  contracting  cones  giving  uniform  velocity  at  the  throat. 
Although  this  design  was  specifically  for  circular  cones,  the  resulting 
square  inlet  cone  fabricated,  together  with  wire  screen  type  flow  straight- 
eners  provided  the  required  uniform  velocity  profile.   Velocity  profiles, 
taken  at  the  matrix  inlet  over  the  entire  flow  range,  verified  required 
uniformity.   The  test  section  was  placed  upstream  of  the  air  metering 
device  so  that  full  advantage  of  the  uniform  velocity  profile  could  be 
realized.   Uniform  temperature  profiles  within  +0.5  °F  were  obtained  by  a 
wire  screen  type  heater  explained  in  detail  in  APPENDIX  B.   The  latter  part 
of  the  first  idealization  pertaining  to  the  thermal  conductivity  of  the 
matrix  cannot  be  achieved  precisely;  however,  the  effects  of  deviating 
from  it  can  be  made  negligible  in  the  solid,  if  resistance  to  heat  transfer 
is  small  within  the  solid  compared  to  the  heat  transfer  resistance  be- 
tween solid  and  gas.   In  other  words,  "transverse  conduction"  effects  have 
been  neglected;  whereas,  the  effects  of  longitudinal  conduction  are  consid- 
ered with  the  use  of  X.  >  the  conduction  parameter.   In  a  matrix  where 
adequate  mixing  takes  place,  no  criterion  is  necessary  for  "transverse 
conduction"  effects. 

The  second  idealization  is  self  explanatory,  i.e.,  the  fluid  is  re- 
stricted to  a  gas. 

Idealization  (3)  for  the  thermal  properties  can  be  guaranteed  if  the 
temperature  differences  are  kept  small.   A  maximum  air  temperature  varia- 
tion from  70  ° F  to  90  ° F  results  in  a  maximum  deviation  from  the  arith- 
metic average  in  the  order  of  +  1.5%  in  viscosity  and  is  negligible  in 
specific  heat. 

The  convective  heat  transfer  coefficient  presented  as  the  final  result 
of  this  experiment  is  an  overall  quantity;  accordingly,  idealization  (4) 
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is  satisfied. 

Satisfying  the  fifth  idealization  of  a  step  change  in  fluid  tempera- 
ture is  perhaps  the  most  rigorous  requirement  of  all.   Several  investiga- 
tors have  attempted  to  meet  this  requirement  by  various  schemes.   The 
technique  of  Bannon  £lj  and  Theoclitus  and  Eckrich  (jEOJ  appears  to  be 
most  satisfactory;  therefore,  it  was  the  one  selected  for  this  investiga- 
tion. 

Figure  1  displays  a  schematic  of  the  experimental  set  up.   Figure  2, 
is  a  photograph  of  the  complete  test  apparatus.   The  system  is  comprised 
of  an  air  source,  flow  straightener  screen  section,  heater  section,  matrix 
test  section,  air  metering  device,  and  pressure  and  temperature  measuring 
systems.   Each  of  these  are  discussed  at  length  in  APPENDIX  B. 

Air  was  drawn  through  the  "wind  tunnel"  by  a  Spencer  Turboblower. 
The  0.0031"  diameter  electric  wire  heaters  upstream  of  the  matrix  were 
energized  so  as  to  heat  the  air  to  roughly  20  degrees  above  ambient.   The 
upstream  and  downstream  fluid  temperatures  at  the  matrix  test  section  were 
monitored  to  ascertain  when  the  matrix  had  been  heated  to  a  uniform  tempera- 
ture.  When  this  condition  was  satisfied,  the  heaters  were  quickly  shut 
down  and  the  time-temperature  history  of  the  entering  and  exit  air  at  the 
matrix  was  recorded  on  a  Minneapolis  Honeywell  "Brown"  Recorder.   The 
recorded  curve  is  the  experimental  cooling  curve  from  which  the  maximum 
slope  was  determined. 

The  mass  rate  of  flow  of  air  was  measured  by  an  ASME  orifice  meter 
with  pressure  taps  one  diameter  upstream  and  one  -  half  diameter  down- 
stream, using  a  range  of  orifice  plates. 

Flow  friction  performance  of  the  matrix  was  obtained  by  measuring  the 
isothermal  pressure  drop  across  the  matrix  and  the  mass  rate  of  flow  of 
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air  through  the  matrix.   Pressure  taps  upstream  of  the  matrix  and  upstream 
of  the  orifice  meter  were  used  to  record  static  pressure  entering  the 
system  and  at  the  orifice  section  respectively. 

Data  reduction  relations  are  discussed  in  APPENDIX  A.  A  sample 


calculation 
The  fol 
(a 
(b 
(c 
(d 
(e 
(f 
(g 
(h 

(i 


or  one  test  run  is  presented  in  APPENDIX  D. 

owing  data  was  recorded  for  each  test  run: 

Orifice  diameter  (d  ) 

o 

Entering  static  pressure  (P  ) 

s 

Pressure  drop  across  the  matrix  (  A  ?m) 

Static  pressure  upstream  of  orifice   (P  ) 

Pressure  drop  across  the  orifice   (A.P  ) 

Orifice  temperature   (t  ) 

Atmospheric  pressure   (P   ) 

atm 

Time  -  temperature  cooling  curve 
Recorder  chart  speed 
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4.    DESCRIPTION  OF  SURFACES 

Surfaces  examined  in  this  experiment  fell  into  one  of  two  major  classi- 
fications -  perforated  nickel  plate  matrices  or  offset  rectangular  fin  of 
stainless  steel.   In  the  former  case,  two  of  the  specimens  examined  - 
the  20  c  skew  matrix  and  a  modified  parallel  plate  matrix,  both  of  160/40TV 
perforated  nickel  -  had  previously  been  examined  by  Bannon  [l]  .   These  parti- 
cular surfaces   were  selected  to  provide  an  additional  check  on  the  system 
reliability. 

The  20  °  skew  matrix  was  formed  by  corrugated  sheets,  which  were  skewed 
relative  to  one  another  to  prevent  meshing  of  adjacent  sheets  at  an  angle 
of  20°  .   This  matrix  configuration  unlike  both  the  modified  parallel 
plate  and  the  offset  rectangular  fin  matrices  does  not  require  splitter 
plates. 

The  "modified  parallel  plate"  configuration  consists  of  formed  sheets 
separated  by  splitter  plates  to  prevent  nesting.   The  assembled  matrix  is 
a  maze  of  rectangular  passages  with  a  value  of  ot   ,  the  aspect  ratio  of  a 
rectangular  flow  passage,  =  7.   The  following  sketch  depicts  the  ideally 
formed  "modified  parallel  plate"  matrix: 


*-=b/c 


0 
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A  modified  parallel  plate  matrix  of  each  type  of  perforated  nickel  was 
constructed  and  tested  to  examine  the  effects  of  varying  shape,  quantity  and 
size  of  perforations  on  this  configuration. 

The  offset  rectangular  fin  surfaces  were  provided  by  the  courtesy  of 
the  AiRsearch  Manufacturing  Company,  Los  Angeles.   The  method  of  identifi- 
cation employed  with  this  configuration  is  best  illustrated  by  the  follow- 
ing example:  [_ 2  J 

16R  -  .153  -  .143  (0)  -  .004  -  SS 
®©    ©     ©  ©     ©  ® 

^  Number  of  fins  per  inch 

(?)   Type  of  fin  flow  cross-section  (R  =  rectangular) 

©  Fin  height,  inches 

@  Uninterrupted  fin  length  in  the  flow  direction,  inches 

©  Type  of  surface  (  0  ■  offset) 

©  Fin  metal  thickness,  £  ,  inches 

©  Fin  material  (SS  =  stainless  steel,  Al  -  Aluminum) 
Two  matrices  of  this  configuration  were  investigated:  -   the  16R  - 
.153  -  .143  (0)  -  .004  -  SS  and  the  20R  -  . 100  -  .125  (0)  -  .004  -  SS. 
The  offset  -  rectangular  fin  matrix  is  comprised  of  offset  fins,  separated 
by  splitter  plates  (0.004"  thick)  to  prevent  nesting  of  the  fins  which 
would  alter  the  flow  characteristics.   The  isometric  sketch  below  best 
identifies  this  plate  fin  arrangement: 
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5.    PRESENTATION  OF  RESULTS 

Heat  transfer  and  flow  friction  data  for  each  matrix  investigated  are 
exhibited  in  both  tabular  and  graphical  form,  TABLES  I  through  X  and 
Figures  9  through  18,  using  the  Colburn  j  modulus,  Fanning  friction  factor, 
and  Reynolds  number.   In  evaluating  the  Reynolds  number,  the  hydraulic 
diameter,  which  is  defined  as  four  times  the  hydraulic  radius,  was  used. 
The  effects  of  entrance,  exit,  and  flow  acceleration  have  been  considered 
in  the  evaluation  of  f.   (See  Appendix  A). 

Geometric  and  physical  properties  of  the  surfaces  are  given  in  Figures 
3  through  8.   Figures  19  through  23  show  heat  transfer  and  friction  compari- 
sons between  various  matrices  investigated.   In  Figures  24  and  25,  "figure 
of  merit"  type  curves  are  presented.  The  j/f  versus  Reynolds  number  thus 
plotted  gives  an  indication  of  the  required  matrix  flow  frontal  area  for 
a  given  pressure  drop. 

Lastly,  in  Figures  26  and  27,  heat  transfer  power  versus  flow  friction 
power  curves  on  a  unit  area  basis,  evaluated  for  fluid  properties  at  stand- 
ard conditions  of  dry  air  at  500  °F  and  one  atmosphere,  are  represented. 
(See  Appendix  A). 
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6.    EXPERIMENTAL  UNCERTAINTIES 

Those  items  of  possible  error  arising  from  lack  of  exact  duplication 
of  the  idealizations,  boundary  and  initial  conditions  are  discussed  in 
the  "Experimental  Technique"  portion  of  this  report  and  no  further  attempt 
was  made  to  evaluate  this  type  of  error. 

The  experimental  errors  which  do  not  fall  into  the  above  area  are 
three  fold  and  are  discussed  below.   These  three  categories  are:   (1) 
uncertainty  of  physical  constants,  (2)  inaccuracy  in  the  geometric  meas- 
urements and  (3)  instrumentation  inaccuracies.   For  the  determination  of 
the  uncertainties  reported,  the  method  of  Kline  and  McClintock  [_12J  was 
used. 

(1)   Values  for  the  physical  constants  required  were  obtained 
from  [4 ^  i  ^5 1  and  \_6j   .   The  uncertainties  in  these  values,  as  best 
as  can  be  determined,  are  listed  below: 


c 
s 

+ 

0.57. 

c 

p 

+ 

0.57. 

NPr 

+ 

2.0% 

k 
s 

+ 

0.5% 

JU 

+ 

1.0% 

It  is  noteworthy  that  the  accuracy  for  the  value  of  c  depends  upon 

the  metal  used.   Where  pure  metals  are  used  the  accuracy  quoted  is  feasible; 

however,  this  is  not  the  case  with  alloys.   Thus  the  tolerance  of  c  can  be 

a  major  cause  for  inaccuracy  ^13j  . 

(2)   Due  to  inconsistencies  in  construction  of  the  matrices, 

errors  in  lineal  dimensions  are  considered  to  be  +  0.5%.   Inasmuch  as  the 

weight  of  the  matrix  can  be  determined  as  accurately  as  desired,  the  error 

in  W  is  considered  negligible.   Based  on  this,  the  uncertainties  in  geo- 
s 

metric  measurements  are  considered  to  be  as  follows: 
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A,  Afr,  Ac,  Ag,  Ak,      +  1.0% 

L       +  0.5% 

W       <  0.1%,  negligible 
s 

(3)   Instrumentation  inaccuracies  were  essentially  those 
obtained  in  pressure  measurement.   With  the  exception  of  the  temperature 
measurement  at  the  orifice,  all  temperatures  recorded  were  temperature 
differences  which  were  recorded  in  inches  (See  Appendix  A).   The  tempera- 
ture measurement  at  the  orifice  was  obtained  by  a  copper-constantan  shield- 
ed thermocouple  which  was  read  by  a  Leeds  and  Northrup  portable  potentio- 
meter.  Assuming  adequate  manufacturer's  calibration  of  the  thermocouple 
wire,  the  estimated  possible  error  is  +  one-half  of  the  smallest  division 
on  the  potentiometer,  or  +  0.0025  millivolts  which  is  approximately  +0.1 
°F. 

Since  the  range  of  pressures  varied,  requiring  several  different  mano- 
meters, the  maximum  uncertainty  encountered  was  used  in  the  uncertainty 
analysis.   These  values  are  listed  below: 

PQ  +  1.25% 

AP  +  .    1.25% 

*— *    o  — 

Z^P  +  1.70% 

m  "~ 

P        +        0.0005"  Hg  (negligible) 
atm 

One  final  parameter  for  which  uncertainty  must  be  determined  before 

an  analysis  of  the  experimental  results  can  be  made  is  N  .   Inasmuch  as 
'  r  Tu 

Colburn  j  factor  is  not  a  linear  function  of  Reynolds  number,  uncertainty 
for  high  and  low  N  values  was  considered.   Using  Figure  30  and  the  un- 

certainty  in  Maximum  Slope  as  2.07.  these  values  were  determinable. 

dNTu 

For  Nm  -  3.0,      X   =  0        Tj =  +  7.2% 

Tu  wTu 
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For  N     -  25.0,    \  =  0.05 

Tu         * 


dN, 


Tu 


=  +   10.07. 


Tu 


Now,  with  the  aid  of  reference  J_12j  ,  the  uncertainties  in  j,  m,  N  ,  and  f 

-iVSl 


were  evaluated  based  on  20:  1  odds. 


^R  " 


where  w  ,  w„,  ....w  ■  uncertainty  interval  based  on  certain  odds, 

1*   2       n  '  * 

and  w  is  the  interval  for  the  result. 
R 

Considering  the  evaluation  of  the  Colburn  j  factor: 


then: 


2/3 


j     =     N, 


Tu        A 


Pr 


Wj    = 


^f&4+  (fc*&( 


'/, 


J  'l 


(^K^V^t(Hi^^^  + 


^^ruAcNPr' 


*C) 


Dividing  by  j  yields: 


Wj    s 


Wntju. 

s.  Nru 


VJi 


+ 


W*c  \    +    / 


3     HPrVj 


-ife 
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For  the  N_       =     3.0,       X    =  0  case: 
Tu  ' 


^   =  [(.o^V(.o,)\(.o,)\(t)  J  -  -0Mt 


W,' 


:S  s  7.5% 
j 


For   the  N         =25.0,       X    =  0.05 
Tu 


case: 


^[(.,)\(.oo\(.ooMi?)^ 


1018 


*i  s  |O.Z% 

j 


By  similar  analysis  the  uncertainties  at  20:1  odds  were  determined 
for  the  following: 


m 


N 


R 


+ 

1.07. 

+ 

2.3% 

+ 

4.37. 
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7.    DISCUSSION  OF  RESULTS 

In  order  to  verify  the  reliability  of  the  data  obtained  by  the  test 
apparatus  used,  several  surfaces  on  which  data  was  currently  available 
were  tested.  Two  of  these  matrices,  the  160/40TV  perforated  nickel 
"parallel  plate"  and  20°  skew,  were  evaluated  by  Bannon  [l]  .   These  two 
matrices  were  originally  evaluated  by  the  same  transient  test  technique 
on  basically  the  same  test  apparatus;  therefore,  satisfactory  comparison 
of  the  data  provided  the  necessary  confidence  in  the  test  rig.   Figures 
9  and  16,  as  well  as  TABLES  I  and  VIII,  present  the  results. 

The  two  offset  rectangular  fin  matrices  tested  were  essentially  tested 
by  Briggs  in  [2  J  with  the  exception  of  slightly  thicker  splitter  material 
and  using  aluminum  in  lieu  of  the  stainless  steel.   The  data  in  £2]  was  ob- 
tained by  a  steady  state  "steam-to-air"  test  technique.   Although  this 
technique  permitted  investigation  aty  much  higher  values  of  Reynolds  number, 
a  sufficient  overlap  in  the  Reynold's  range  existed  so  that  a  comparison 
was  possible.   A  good  comparison  for  the  20R  -.100  -   125  -  .004  -  SS 
surface  for  Colburn  j  factor  exists;  however,  the  friction  factor  re- 
ported in  ^2  J  differs  noticeably.   Results  are  shown  in  Figure  17  and 
TABLE  IX.   For  the  16R  -  . 153  -  .143  -  .004  -SS  surface,  the  values  for 
Colburn  j  factor  for  Reynolds  number  <  1000  are  in  excellent  agreement 
with  those  presented  in  ^2J  .   In  the  transition  range,  the  j  values  re- 
ported herein  are  noticeably  different;  however,  the  values  for  friction 
factor  agree  nicely  throughout  the  entire  duplicated  Reynolds  number  range. 
It  is  important  to  note  here  that  by  the  "maximum  slope"  method  of  analysis 
an  inherently  large  error  is  possible  in  the  low  N   region  (1.5  f=  N„,  S: 

3.5)  r 71  .   The  size  of  this  error  in  N   relative  to  an  error  in  maximum 
L  J  Tu 

slope  can  best  be  appreciated  from  Figure  30.   Consider,  for  example,  a 
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2%  error  in  maximum  slope  with  an  N   value  of  2.3  and  X.  ~  0.   The 

Tu 

percent  error  in  N  ,  from  Figure  30,  is  of  the  order  of  207..   For  values 
of  N   between  2.0  and  2.3,  the  percent  error  has  not  been  evaluated.   The 
1.5  ^  NTu  ^  3.5  region  coincides  with  the  Reynolds  range  for  which  the 
reported  j  values  differ  from  [l  ]  ,  and  may  account  for  the  difference 
noted.   The  results  for  this  matrix  appear  in  Figure  18  and  TABLE  X. 

In  spite  of  the  areas  of  disagreement,  the  correlation  was  such  that 
the  data  for  the  j  in  the  Reynolds  range  < 1000  and  the  friction  factor 
data  for  these  matrices  are  considered  reliable.   It  is  noteworthy  that 
the  surface  tested  in  [2]  is  specified  as  15.76R  -  . 153  -  .143  -.004  -  Al, 
so  that  the  matrix  evaluated  herein  may  only  be  an  approximate  duplicate. 
A  comparison  of  the  relative  performance  for  these  matrices  is  presented 
in  Figure  23.   Additionally,  Figure  25  illustrates  j/f  vs.  NR,  a  flow 
area  "goodness"  factor. 

The  perforated  nickel  matrices  investigated  covered  a  wide  range  of 
percent  open  area,  12  to  50%.   Both  slotted  and  round  perforated  nickel 
surfaces  were  evaluated  for  the  "modified  parallel  plate"  geometry.   The 
perforated  nickel  geometric  properties  are  given  in  APPENDIX  C  for  the 
various  types  investigated.   In  order  to  make  a  reasonable  comparison  of 
the  performance  among  these  matrices,  the  same  frontal  area  and  flow  length 
were  specified.   Results  are  presented  in  Figures  9  and  10  and  listed  in 
TABLES  I  and  II  for  the  slotted  perforated  nickel  matrices;  those  for  the 
round  perforated  matrices  are  given  in  Figures  11,  12,  and  13  and  tabu- 
lated in  TABLES  III,  IV,  and  V. 

In  order  to  establish  a  meaningful  fiducial  point,  a  solid  nickel 
matrix  of  the  identical  material  (obtained  from  the  same  manufacturer) 
with  the  same  frontal  area  and  flow  length  as  the  perforated  matrices  was 
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fabricated.   The  performance  of  this  matrix  is  illustrated  in  Figure  15 
with  tabular  values  given  in  TABLE  VII. 

The  best  performing  slotted  perforated  nickel  matrix  was  the  160/40Q. 
The  50G  round  perforated  nickel  (507.  open  area)  was  the  best  surface  in- 
vestigated among  the  "parallel  plate"  nickel  matrices. 

As  another  approach  to  evaluating  the  Colburn  j  factor,  J'  was  con- 
sidered.  The  parameter  j'  is  based  on  the  plane  surface  area,  A*,  as  if 
no  perforations  were  present;  whereas, j  is  based  on  the  heat  transfer  sur- 
face area,  A,  which  includes  the  solidity  correction,  s.   The  value  of  j' 
is  readily  attainable  from  the  relationship 

j'  =  j  x  s     ,         where  s  =  A/ A*.   (See  Appendix  D) 
The  determination  of  the  solidity  is  fully  described  in  APPENDIX  C.   Figure 
13  for  the  50G  perforated  nickel  matrix  includes  the  values  for  j'  versus 
Reynold's  number  in  addition  to  the  f  and  j  values. 

A  comparison  of  the  performance  of  the  two  slotted  perforated  matrices 
together  with  the  solid  nickel  matrix  is  shown  in  Figure  19.   A  similar 
comparison  for  the  round  perforated  nickel  matrices  is  given  in  Figure  20. 

It  is  readily  apparent  that  the  perforated  material  yields  a  higher 
heat  transfer  characteristic.   Inasmuch  as  the  friction  factor  also  shows 
an  increase,  a  more  appropriate  comparison  can  be  made  from  the  flow  area 
"goodness"  factor  curve  in  Figure  24.   The  heat  transfer  power  (h   )  versus 

5  LU 

flow  friction  power  (ECTn)  per  unit  area  given  in  Figures  26  and  27,  for  the 
slotted  and  round  perforated  nickel  matrices  respectively,  is  another  "good- 
ness" factor  evaluation.   Here  the  information  for  each  of  the  matrices  is 
presented  based  on  standard  conditions  for  fluid  properties  and  a  common 
hydraulic  diameter.   (See  APPENDIX  A). 
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So  that  an  appreciation  for  the  perforated  material  might  be  realized 
in  direct  transfer  type  of  heat  exchanger  application  in  addition  to  its 
use  in  matrix  type  for  rotary  regenerators,  a  "parallel  plate"  matrix  with 
the  160/40TV  perforated  nickel  fins  and  solid  nickel  splitter  plates  was 
fabricated.   The  performance  of  this  matrix  is  given  in  graphical  and 
tabular  form  in  Figure  14  and  TABLE  VI.   A  comparison  with  the  solid 
nickel  matrix  is  made  in  Figure  21.   Figure  24  contains  the  j/f  presenta- 
tion and  Figure  26  shows  the  h vs.  E  for  this  matrix.   Here  again  is 

STD      STD  6 

a  noticeable  increase  in  heat  transfer  performance,  but  with  a  negligible 
increase  in  friction  factor. 

For  the  160/40TV  perforated  nickel  material,  a  20°  skew  matrix  was 
compared  in  performance  against  the  "modified  parallel  plate"  to  see  how 
these  two  configurations  compare  with  one  another.   The  separate  perform- 
ance for  each  has  been  mentioned  previously;  however,  the  direct  comparison 
is  shown  in  Figure  22.   As  with  the  other  nickel  surfaces,  the  flow  area 
"goodness"  factor  is  represented  in  Figure  24. 

In  the  results  of  this  investigation  are  several  items  which  are  general 

in  nature,  but  nonetheless  extremely  important.   In  the  laminar  Reynolds 

number  range,  friction  factor  is  inversely  proportional  to  Reynolds  number. 

By  Reynolds  analogy,  it  was  anticipated  that  the  Colburn  j  factor  should 

likewise  be  inversely  proportional  to  Reynolds  number  in  this  region. 

However,  the  experimentally  obtained  values  of  j  in  the  low  laminar  Reynolds 

range  utilizing  the  equivalent  conduction  length  (defined  in  APPENDIX  C) 

together  with  TABLE  XI  or  Figure  29  yielded  lower  values  for  j.   For  this 

reason,  a  re-evaluation  of  the  conduction  parameter,  N   ,  and  maximum  slope 

Tu 

data  points  for  Figure  29  was  performed.   The  equivalent  conduction  para- 
meter X  ,  ,  was  computed  and  the  N   value  was  determined  by  utilization  of 
k  Tu 

the  extrapolated  j  vs.  N  curve,  since  j  is  proportional  to  N   .   This  new 
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determined  value  of  N   was  considered  to  be  the  actual  N   for  the  parti- 

Tu  Tu 

cular  "maximum  slope"  and  conduction  parameter.   Values  determined  in  this 
fashion  are  asterisked  in  TABLES  I  through  VI  and  VIII.   The  tabulated 
conduction  parameter  values  in  these  tables  are  the  equivalent  conduction 
parameter  (See  APPENDIX  D) . 

The  second  item  is  that  of  an  equivalent  conduction  length  which  was 
quoted  earlier.   In  the  perforated  fin  material  the  slots  or  round  holes, 
as  the  case  may  be,  provide  for  a  winding  heat  flow  path.   It  was  appreci- 
ated that  the  matrix  flow  length,  L,  was  not  the  true  conduction  path  but 
somewhat  less  than  actual.   A  mean  conduction  path  length  evaluation  is 
presented  in  APPENDIX  C.   The  conduction  parameter,  A  ,  based  on  the  flow 

length,  L,  was  evaluated,  and  by  simply  multiplying  by  the  ratio  L/L   (also 

k 

given  in  APPENDIX  C),  the  value  of  A   was  determined. 

k 

A  third  item  is  the  inherent  error  in  N   in  the  range  1.5—  N„  —  3.5 

Tu  Tu 

which  was  mentioned  earlier.   Howard  [7  J  published  an  error  curve  for  N   , 

Tu 

including  slope  vs.  N  ,  see  Figure  30. 

Tu 

Fourthly,  the  entrance  and  exit  length  effects  for  a  matrix  with  a 
small  hydraulic  radius  are  small.   In  the  evaluation  for  Fanning  friction 
factor  (APPENDIX  A),  the  effects  of  exit  and  entrance  length  were  consider- 
ed using  values  of  K  and  Ke  from  Figures  37,  38,  or  39,  as  appropriate. 
For  the  perforated  surfaces,  the  K  values  were  taken  corresponding  to 
laminar  flow  conditions;  however,  in  the  offset  rectangular  fin  geometry 

the  values  were  taken  from  those  Figures  for  N  =  00  in  accordance  with 

K 

reference  {llj    •   The  offset  fins  never  have  velocity  profiles  that  are 
fully  established,  that  being  a  purpose  in  the  design.   The  same  calcula- 
tions were  repeated  on  the  offset  rectangular  fins  using  the  K  values  for 
laminar  flow.   The  friction  curves  plotted  for  the  offset  rectangular  fins 
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were  coincident  with  those  obtained  earlier;  however,  some  change  was  noted 

in  tabulated  values. 

A  final  item  worth  mentioning  is  again  concerned  with  the  low  Reynolds 

range.  A  slight  alteration  in  the  conduction  parameter  value  has  a  marked 

effect  on  the  N   value  in  this  range,  which  is  apparent  from  Figure  29. 
Tu 

For  this  reason,  when  utilizing  a  nickel  material,  or  any  other  material 

with  a  wide  range  of  specified  values  for  thermal  conductivity,  it  is 

mandatory  that  the  appropriate  value  of  k  be  obtained  for  calculating 

s 

purposes.   The  value  of  \.  is  directly  proportional  to  this  k  value, 

s 


i.e. , 


A.  =k- A 


_s S 

L  m  cp 
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8.    CONCLUSIONS 

Reliable  data  is  obtainable  from  the  experimental  apparatus  based  on 
a  comparison  with  those  surfaces  previously  evaluated  elsewhere. 

A  very  definite  improvement  in  performance  is  achieved  with  the  utiliza- 
tion of  a  perforated  material.  Kays,  in  reference  [_l6]    ,  postulated  the 
improved  heat  transfer  characteristic  without  a  "large  amount  of  form  drag 
so  characteristic  of  high  performance  surfaces".   The  perforations  apparent- 
ly disturb  the  heat  transfer  sub  boundary  layer  so  that  tubulence  increases 
and  the  resulting  high  convection  conductance  is  realized.   Inasmuch  as  the 
perforations  are  fine,  the  downstream  edge  is  considered  to  be  in  the  wake 
of  the  separated  region.   The  perforated  nickel  matrices  investigated  here- 
in confirm  Rays  hypothesis.   Unfortunately,  inasmuch  as  the  perforated  sur- 
face investigated  in  [loj  is  greatly  different  in  geometric  and  physical 
properties,  no  meaningful  comparison  with  those  tested  in  this  investiga- 
tion could  be  made. 

As  stated  previously,  the  experimentally  obtained  Colburn  j  values  in 

the  low  laminar  Reynolds  range  using  TABLE  XI  or  Figure  29  were  lower  than 

j  values  obtained  by  extrapolating  the  j  versus  N  curve.   Inasmuch  as  fric- 

R 

tion  factor  and  Colburn  j  factor  are  both  inversely  proportional  to  Reynolds 
number  in  the  laminar  region  owing  to  Reynolds  analogy,  the  values  for 
Colburn  j  factor  extrapolated  from  the  experimental  curves  are  considered 
more  reliable  in  the  low  laminar  range. 

Effects  of  longitudinal  conduction  have  been  considered  in  the  analysis 
of  data.   For  the  perforated  material  matrices,  corrections  for  the  actual 
conduction  path  length  and  for  the  area  removed  by  the  perforations  have 
been  included. 

The  best  performance  for  a  perforated  "parallel  plate"  matrix  was  at- 
tained by  the  50G  type  material.   For  the  offset  rectangular  fin  surfaces, 
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the  20R  -  .100  -  •125  -  -004  -  SS  matrix  was  superior, 
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9.    RECOMMENDATIONS  FOR  FURTHER  STUDY 

1.  Due  to  the  inherent  error  in  the  high  Reynolds  number  range  where 

1.5  ^  N   ^3.5,  a  more  reliable  means  for  determining  N   values  should 

be  investigated.   Perhaps  pursuit  of  Hausen's  cyclic  operation  [_9j  ,  a 

steady  state  periodic-flow  technique,  is  appropriate  in  that  this  method 

measures  effectiveness  from  which  N   can  be  determined  subsequently. 

Tu 

Since  the  results  obtained  using  Schumann's  analysis,  even  with  longitudin- 
al conduction  effects  considered,  at  the  low  Reynold's  range  appears  in- 
adequate; the  use  of  another  testing  method  may  prove  fruitful. 

2.  Inasmuch  as  the  best  performing  surface  for  perforated  nickel 
was  the  50G  type  (50%  open  area),  other  types,  such  as  the  40G  which  has 
a  slightly  larger  hole  diameter,  with  same  plate  thickness  and  percent 
open  area  as  the  50G  may  prove  worthy  of  test. 
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Matri      rial 

Specific  Heat  (c3)    •     F 

Thennal  Conductivity  (ks    -        f 

Material  Thickness,  inches 

Total  Heat  Transfer      ,A)  ft^ 

Fn  : 

Total.  Conduction  Area  Uk)   ft 

Free  Flow  Area  Uc) 
; rix  Volume  (Vm)   :  ' 
-be  Den^i    f  :n)  lb/f' 

Hydraulic  Diameter  I 

Compactness  (@) 

Porosity  (p) 
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160, .. 
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160/40   v 
Perf  .Mck 
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16.2667 
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0.011588 

63.6 

59.2 
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Matrix  Material 

Specific  Heat  (cs)  Btu/lb  °F 

Thermal  Conductivity  (ks)  Btu/  hr  ft  °F 

Material  Thickness,  inches 

Total  Heat  Transfer  Area  (A)   ft2 

Frontal  Area  (Afr)   ft2 

Total  Conduction  Area  (A  )  ft 

Free  Flow  Area  (Ac)  ft2 

Matrix  Volume  (Vm)   ft3 

Matrix  Density  (  ^  m)  lb/ft3 

Hydraulic  Diameter  (DH) 

Compactness  ( 0 )  ft2/ft^ 

Porosity  (p) 


Solid  Nickel 
0.1065 

38.7 
0.0020 
20.1875 
0.06953 

o.oiou 

.     939 
0.011538 
77.7 

0.001961 
1742.1 
0.854 


Fif^ure  6.  Geometric  and  Physical  P. 

Solid  Nickel  "  ix 
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Matrix  Material 

Specific  Heat  (cs)  Btu/lb  F 

Thermal  Conductivity  (ks)  Btu/  hr  ft     F 

Material  Thickness,  inches 

Total  Heat  Transfer  Area  (a)  ft2 

Frontal  Area  (Afr)  ft2 

Total  Conduction  Area  (A^)  ft2 

free  Flow  Area  (Ac)  ft2 

Matrix  Volume  (Vm)  ft3 

Matrix  Density  (^m)  lb/ft3 

Hydraulic  Diameter  (D^)  ft 

Compactness  ( f )  ft2/ft3 

Porosity  (p) 


20  Degree  Skew 
160/40  TV 
Perf .Nickel 

0.1065 

38„7 

0.0022 

15.41 

0.0734 

0.00349 

0.06412 

0.01211 

50.1 

0.00252 

1272.8 

0.873 


Figure  7.  Geometric  and  Physical  Properties  of  a  Slotted 
Perforated  Nickel  (160/40  TV)  20°  Skew  Matrix 


Matrix  Material 

Specific  Heat  (cs)  Btu/lb  F 

Thermal  Conductivity  (ks)  Btu/  hr  ft  F 

Material  Thickness,  inches 

Total  Heat  Transfer  Area  (A)   ft2 

Frontal  Area  (Afr)  ft2 

Total  Conduction  Area  (As)  ft2 

Free  Flow  Area  (Ac)   ft2 

Matrix  Volume  (Vm)   ft3 

Matrix  Density  (  ^  m)  lb/ft 3 

Hydraulic  Diameter  (D^)   ft 

Compactness  (Q)      ft2/ft3 

Porosity  (p) 


2Urt  -.100  - 

16R  -.153  - 

.125-.004-SS 

.143-.C04-SS 

0.120 

0.120 

9.3 

9.3 

0.004 

0.004 

11.625 

9.0625 

0.07083 

0.C7083 

0.01549 

0.01201 

0.05535 

0.05882 

0.0177 

0.0177 

66.3 

50.4 

0.00476 

0.00649 

656.5 

511.8 

0.781 

0.830 

Figure  8.  Geometric  and  Physical  Properties  of  Stainless 

Steel  Offset  Rectangular  Fin  (20R-.100-.125-. 004 
-SS  and  16R-.153-.143-.004  -SS)  Matrices 
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Figure  9»  '  Surface  Heat  Transfer  and  Friction  Data 

"Parallel  Plate",  160/40  TV  Perforated  Nickel 
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Figure  10.  Surface  Heat  Transfer  and  Friction  Data 

"Parallel  Plate",  160/40  Q  Perforated  Nickel 
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Figure  11.  Surface  Heat  Transfer  and  Friction  Data 

"Parallel  Plate",  125  M  Perforated  Nickel 


Figure  12.  Surface  Heat  Transfer  and  Friction  Data 

"Parallel  Plate",  125  P  Perforated  Nickel 
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Figure  13.  Surface  Heat  Transfer  and  Friction  Data 
"Parallel  Plate",  50  G>  Perforated  Nickel 
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Figure  14.  Surface  Heat  Transfer  and  Friction  Data 

"Parallel  Plate",  160/40  TV  Fine  with  Solid 
Nickel  Splitter  Plates 


Figure  15 •  Surface  Heat  Transfer  and  Friction  Data 
"Parallel  Plate",  0.002"  Solid  Nickel 
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Figure  16.  Surface  Heat  Transfer  and  Friction  Data 
20*  Skew,  160/40  TV  Perforated  Nickel 
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Figure  17.  Surface  Heat  Transfer  and  Friction  Data 

Offset  Rectangular  Fin,  20R  -.100  -.125  -.004  -  SS 


Figure  18.  Surface  Heat  Transfer  and  Friction  Data 

Offset  Rectangular  Fin,  16R  -.153  -.H3  -.004 


-  SS 
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Figure  19,  Comparison  of  Slotted  Perforated  vs.  Solid 
Nickel  "Parallel  Plate "  Matrices 
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Figure  20.  Comparison  of  Round  Perforated  vs.  Solid 
Nickel  "Parallel  Plate"  Matrices 
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Figure  21.  Comparison  of  Slotted  Perforated  Fin  with  Solid 

Nickel  Splitters  vs.  Solid  Nickel  "Parallel  Plate" 
Matrices 
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Figure  22.  Comparison  of  Slotted  Perforated  Nickel  "Parallel 
Plate"  vs.  Slotted  Perforated  Nickel  20°  Skew 
Matrix 
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Figure  23.  Comparison  of  Offset  Rectangular  Fin  Matrices 


51 


TO 
<D 
O 

u 


JS 

-p 

O 

<M 

TO 
U 
O 

-P 

3 

TO 
TO 
Q 


O 
O 
O 


2 

<: 

8 

si 


:: 


Figure  25.  Flow  Area  Goodness  Factors  for  the  Offset 
Rectangular  Fin  Surfaces 
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FIGURE  28.  GENERALIZED     COOLING     CURVE 
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Figure  31.  Matrix  Holder 
(Upstream  View) 
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Figure  34.     Ni chrome  Wire  Heater 
(Upstre?im  View) 


Figure  37.  Entrance  and  Exit  Pressure  Loss  Coefficients 
for  Triangular  Cross  Section  Matrix  with 
Abrupt-Contraction  Entrance  and  Abrupt- 
Expansion  Exit 
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Figure  38 •  Entrance  and  Exit  Pressure  Loss  Coefficients 
for  Parallel  Plate  Cross  Section  Matrix,  with 
Abrupt-Contraction  Entrance  and  Abrupt- 
Expansion  Exit 
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Figure  39*  Entrance  and  Exit  Pressure  Loss  Coefficients 

for  Multiple  -  Rectangular  Cross  Section  Matrix, 
with  Abrupt-Contraction  Entrance  and  Abrupt- 
Expansion  Exit 
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APPENDIX  A 
DATA  REDUCTION  RELATIONSHIPS 
Summarized  herein  are  those  data  reduction  relationships  which  are 
of  prime  importance  in  obtaining  the  heat  transfer  and  flow  friction 
characteristics  as  presented  from  the  raw  data  taken.   Before  the  data 
on  the  matrices  tested  can  be  reduced,  or  a  comparison  of  the  different 
matrices  can  be  made,  it  is  essential  to  obtain  accurate  geometrical  dimen- 
sions.   The  geometrical  factors  of  interest  are  the  porosity,  p,  the  hy- 
draulic diameter,  D  ,  and  the  area  compactness,  (3  .   Determination  of 
any  of  these  fixes  the  third  since  it  may  be  shown  that 


(3  =  ±E 


[14] 


(A-l) 


The  definitions  of  these  terms  are: 


(1) 


Free  Flow  Area 
p  ■  Frontal  Area 


(A-2) 


4  x  Free  Flow  Area 

(2)      D  ■  Heat  Transfer  Surface  Area 
n 


Heat  Transfer  Surface  Area 
(3)      @  =  Frontal  Area  x  Flow  Length 


4AcL 
A 

f\_ 

A^rL 


(A-3) 


(A-4) 


MAXIMUM  SLOPE 

As  presented  earlier,  the  maximum  slope  of  the  generalized  cooling 

curve  is  a  unique  function  of  N_  and  /L  . 
^  Tu 


a 


=  p   (^ru.A) 


(A-5) 


J  rna->c 
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T*  and  NT  have  both  been  previously  defined,  but  for  convenience 
the  relationships  are  restated: 


Wscs 

and 

Ntu 

=  hA 

rnCp 

Therefore: 

/     •                    y 

and 


<*(^/ntu.)  =  *cf/Vs)  de  =  (c/cs)  de 


(A-6) 


Furthermore: 


+    1 


[7] 


and 


t^-ti 


1 


■tu-tij      -tfl-tv 


[t«l-tft] 


(A- 7) 


Combining  (A-6)  and  (A- 7) 


[ 


'^-ti 


[r/NTu] 


c 


■fcft-ti      de 


max 


(A-8) 


From  Figure   28: 


(1) 


(2)       2L 


iaf-t^-  tfi1 
de 


X 


max 


Chart   speed 


=  ae   »  s 
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ec 


(3)  y  =  <*("t£z  -  "t.fl  ),  Inches 

(4)  M  =  (-tfl  -  -t'x  ),   inches 
Combining  these  determined  values  into  (A-8)  yields: 


-   Cg    1   y  .  Chart  speed   (A-9) 

C   '  M  '  * 


max 


With  this  value  of  maximum  slope  and  /\_  =  — 


ksAs 


rnCpL 


,  enter  TABLE 


XI  or  Figure  29  to  obtain  a  corresponding  N   value. 

Tu 


MASS  RATE  OF  FLOW; 

The  mass  rate  of  flow  is  calculated  from  the  method  given  in  reference 
|_16j  with  necessary  rearrangements  in  order  to  reduce  unnecessary  hand 
calculations  and  to  gain  the  full  benefits  of  Murdock  Il5 J  : 


m  =  359Kd, 


Fa  Y 


{m\T 


(A-10) 


where:    K  ■  flow  coefficient,  including  velocity  of 

c 

approach  ■  .<    <r>4=^~   with  C  ■  orifice 

coefficient  of  discharge  (obtained  from  [_15J  ) 

and  Q    =   ratio  of  orifice  diameter  to  internal 

pipe  diameter,  d  /d 
o 

d  =  orifice  diameter,  inches 
o 

F  =  thermal  expansion  factor 
a 


Y  =  expansion  factor 

Y=  specific  weight  of  fluid  flowing  =  ~ 

with   P  in  #/ft 


P_ 
RT 
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R  -  53.35   ft   -  #/lb  °R  for  air 

T  -  t(°F)  +  459.7  -  degrees  Rankine  (°R) 

APQ  -  Pressure  drop  across  the  orifice,  inches  H.O 

substituting  P_  for  V  ,  since  the  magnitude  of  the  local  acceleration  of 

R  T 
gravity,  g,  is  taken  as  being  equal  to  the  magnitude  of  standard  accelera- 

c 

tion  of  gravity,  gc;  and  ^r    _  ~ a,  v  for  K. 
gives: 


Inserting  the  numerical  value  for  R  and  converting  P  to  units  of 


#/in2  yields: 


2. 


A  -  ^9.81C^0  F,Y  "J  APo  .  ^       (a_i2) 

From  Figure  38,  reference  [l6~]    :   F  ^1.00 
From  Figure  40b,  reference  [l 6]  :  10*0.9985^  1.0 
The  pressure,  P,  in  (A-12),  is  hereafter  referred  to  as  P«  and  is 
determined  by  the  following  relationship: 

P.  =  (P  ,   -  P  /  13.6)  0.4912,  #/in2  (A-13) 

j     atm    o 

Where  P     is  the  atmospheric  pressure  in  inches  of  mercury, 
atm 

P  is  the  static  pressure  upstream  of  the  orifice   in 
o 

inches  of  H_0 

13.6  inches  of  H  0/  inches  of  mercury 
2 
2 
0.4912  #/in  /  inches  of  mercury 
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The  temperature,  T,  in  (A-12)  is  the  air  temperature  at  the  orifice 


in  degrees  Rankine  (T  ) : 


TQ  =  tQ  +  459.7  (A- 14) 


Replacing  F  and  Y  in  (A-12)  with  their  numerical  values,  the  equa- 
tion  reduces  to  the  working  form: 


5-89.81  C  d*   I  AP  .  Vz 


■  -   v/  :j^  ak°  *  -^-  (A-i5) 


REYNOLDS  NUMBER 


Reynolds  number  is  defined  as: 


NR  -   ^  (A-16) 


where  G  is  the  mass  flow  velocity  based  on  the  free  flow  area,  A  . 

c 


&- 


no        no 


Ar       P'V  (A-17) 


Substituting  (A-17)  into  (A-16): 


*     /^A^    P  (a-18) 


But  from  (A-l)   VP  =  P   '  therefore:   NR  =  -^_/_i_\   (A-19) 

Equation  (A-19)  shows  that  for  a  given  matrix  and  mass  flow  rate  the 
Reynolds  number  is  inversely  proportional  to  the  compactness.   Substitut- 
ing (A-4)  into  (A-19)  yields  the  working  equation: 

AA  (A.  20) 
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FANNING  FRICTION  FACTOR: 


The  following  sketch  together  with  equation  (A- 21)  describes  the  flow 
system  under  consideration:  1  111 


Flow 


rjwwwwwwwji 


i^wwwwwwwu 

CA\\\\\\\\\\\\\N 

ijwwwwwwwHi 
k\\\\\\\\\\\\\\kl 

i^wwwwwwy 


5=> 


AP     _      G 

Pi 


vl 


2gc       Pi 


(Kc   -  1  -    pZ)    -   2 


*■     entrance 
effect 


flow 
acceleration 


+  f  A, 


'm 


vl 


core 
friction 


-  (1  -  P2  -  Ke) 


exit 
effect 


vl 


(A-21) 


It  is  important  to  note  that  v  =*  v^  and  v,^  v  in  the  above  sketch 
as  the  pressure  changes  (for  gas  flow  heat  exchanger  application)  from 
section  1  to  a  and  b  to  2,  respectively,  are  very  small  relative  to  the 
total  pressure.   Inasmuch  as  the  testing  is  performed  with  air  at  moder- 
ate temperatures  and  pressures,  the  perfect  gas  law  is  applicable  (P=o^RT), 


and  recalling  that  -JL  =  ts.   or   PH    =  -^ 
£  4L     A       L      A 

from  equation  (A-3),  equation  (A-21)  solved  for  f  for  the  isothermal 


Substituting  v  =  -1 


case  becomes: 


\ 


2gc  i 


(  AP 


V  G2  / 


?1+  ?2 


Pi 


+ 


Ke_]    _ 
PjJ 


con't   on 
next  page 
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(con't.) 


(H^U-^K1^2)]^ 


(A-22) 


P  ■+-  P 
where   ^  m  =  -±-L— — ~.      with  the  upstream  values  subscripted  by  1 

and  the  downstream  subscripted  by  2.   K  is  the  entrance  coefficient  and 

Ke  the  exit  coefficient,  and  both  are  dependent  on  porosity,  shape  of 

the  flow  cross  section  and  the  matrix  Reynolds  number.   Values  of  K  and 

K  are  obtained  from  Figures  37,  38,  or  39,  which  were  reprinted  from 

reference  \_11  J  . 

If  we  consider  an  order  of  magnitude  approximation,  the  first  term  in 

(A-22)  is  by  far  the  greatest  contributor  to  friction  factor  for  small 

pressure  differentials. 


The  approximation,  therefore,  that  ^     -p 


l_=Vm<*Vx*>\ 


reduces  (A-22)  to  the  following 


f  -  [2gc  em(£? )  -  (Kc+Ke)+^  (i  +  P*) 


rH 


(A-23) 


Reference  fl4]  points  out  that  the  relationships  given  here  can  be 
recombined  in  such  a  manner  so  as  to  show  that  f  is  proportional  to  the 
porosity  cubed  and  inversely  proportional  to  the  compactness  (f  oC  p  /  (3  ) 

COLBURN  J  FACTOR: 

Colburn  j  factor  is  defined  as: 

■J-    NstlV2/3    -Jl-     NprV3  (A- 24) 
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Substituting  (A-17)  for  G  and  multiplying  by  A/A  yields 


mCp       ^ 

but  U 


j   -   _LA_     J±     NprV3  (A-25) 


mcp 


=      Wtu 


therefore: 


j  -         NTu  4s-  Npr^3  (A-26) 

A 


Combining  equations  (A-l)  and  (A-3): 

Ac  _  ? 
A   @L_ 


thus: 


j  =  ^ru.NPrV3  JL  X  (a-27) 


L  @ 


It  is  apparent  from  equation  (A-27)  that  the  Colburn  j  factor  is 
directly  proportional  to  porosity  and  inversely  proportional  to  compact- 
ness  (jcc  P/^). 

HEAT  TRANSFER  POWER  AND  FLOW  FRICTION  POWER: 

An  evaluation  of  the  heat  transfer  power  versus  flow  friction  power 
is  of  interest  in  that  it  is  a  measure  of  relative  performance. 

The  heat  transfer  power  per  unit  area  per  degree  temperature  differ- 
ence is   111: 

h  -   Cp£_   J_  (A-28) 

NpT7*     Eh  J   * 
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Equation  (A-28)  with  the  properties  evaluated  at  standard  conditions 
of  dry  air  at  500  F  and  one  atmosphere  becomes: 

hSTD  =  0'02195  (\  )(NR  j  )   Btu/hr  ft2  °F  (A-29) 

where     c  =  0.2477  Btu/lb  °F 
/*-  =  0.0678  lb/hr  ft 
^  =  0.0413  lb/ft3 

Nn   =  0.671 
Pr 

The  flow  friction  power  per  unit  area  is  fill  : 

■&  f4-f  f  ^r3  (A"30> 


2gc   ^  VCH 

When  equation  (A-30)  is  evaluated  at  standard  conditions   (  la    and  o 
values  given  above),  the  flow  friction  power  per  unit  area  at  standard 
conditions  is: 

E._n  =  1.11  x  10-7  (^r-  )  f  (3lSL-\  HP/ft2 

STD  v  Uh7    VlOOO/ 

For  comparison  purposes  the  surface  geometries  were  reduced  to  a 

common  hydraulic  diameter  of  D  =  2  x  10" J  ft.   Figures  26  and  27  show  such 

n 

comparisons  for  the  nickel  surfaces. 
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APPENDIX  B 
DESCRIPTION  OF  EQUIPMENT 

Heat  transfer  data  was  obtained  utilizing  the  "single  blow"  technique 
which  consists  briefly  of  monitoring  the  fluid  temperature  response  at 
the  test  matrix  exit  while  subjecting  the  inlet  to  the  test  matrix  to  a 
step  change  in  fluid  temperature. 

Friction  factor  was  determined  from  pressure  drop  data  obtained  from 
static  pressure  taps  located  in  the  test  section  -  one  at  inlet,  the  other 
at  the  exit. 

Necessary  equipment  to  perform  such  an  experiment  falls  into  one  of  the 
following  categories: 

(1)  Fluid  source 

(2)  Flow  metering  system 

(3)  Temperature  measuring  system 

(4)  Pressure  measuring  system 

(5)  Fluid  heater  system 

(6)  Matrix  holder  and  test  section  casing 

FLUID  SOURCE: 

Air,  the  working  fluid,  was  provided  to  the  test  apparatus  by  connect- 
ing the  rig  to  the  inlet  of  a  30HP,  multistage,  Spencer  Turbo-Compressor 
rated  for  550cfm  operating  on  a  220  volt  a.c.  power  supply. 

FLOW  METERING  SYSTEM:   (See  Figure  1) 

The  flow  metering  system  consisted  of  an  ASME  standard  orifice  section 
constructed  for  d  and  d/2  pressure  taps  in  a  3.08"  inside  diameter  metal  tube 
[l6l.  A  wide  flow  range  was  obtained  by  utilizing  thin  plate  concentric 
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orifices  of  throat  diameter  2.310,  1.971,  1.232,  1.540,  0.462,  and  0.308 
inches  respectively. 

Control  of  the  pressure  drop  across  the  orifice  was  maintained  by  a 
gate  valve  downstream  of  the  apparatus,  ahead  of  the  compressor  inlet  and 
another  gate  valve  on  the  compressor  inlet.   The  blast  gate  on  the  compressor 
discharge  was  pre-set  so  that  the  unit  could  not  be  operated  beyond  the  full 
load  rating  of  the  turbo. 

This  turbo-compressor  had  many  distinct  advantages  over  the  one  used  in 
reference  ^lj  :   some  of  the  more  pertinent  being  (1)  a  550  cfm  capacity; 

(2)  essentially  constant  pressure  throughout  the  range  of  the  machine; 

(3)  the  load  decreases  as  less  air  is  used  so  that  the  problem  of  stalling 
at  low  flow  rates  is  non-existent. 

TEMPERATURE  MEASURING  SYSTEM;   (See  Figure  1) 

The  majority  of  the  temperatures  were  measured  with  iron-constantan 

thermocouples*   These  thermocouples  use  iron  for  the  positive  conductor 

and  constantan  for  the  negative  conductor.   Measuring  T  is  a  group  of  five 

thermocouples  in  series,  bound  together  and  insulated  from  each  other  by 

teflon  tape  and  inserted  in  a  1/8"  diameter  aluminum  tube  to  serve  as  a 

radiation  shield.   A  small  aperature  was  cut  in  the  inserted  tube  facing 

upstream.   T  and  1,    are  five  wire  thermocouple  grids,  connected  in  series, 

mounted  permanently  in  the  test  section  casing.   T  is  a  movable  five  wire 

3 

series  connected  thermocouple  grid  located  in  the  matrix  holder  so  that  it 
is  adjacent  to  the  downstream  side  of  the  matrix.   The  purpose  of  connecting 
the  thermocouples  in  series  was  to  magnify  the  emf  output  so  that  the  instru- 
ment sensitivity  to  small  temperature  changes  would  be  enhanced.   With  T_ 

bucked  against  T  ,  it  is  possible  to  determine  the  temperature  uniformity 
4 
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across  the  matrix.   T,  versus  T  indicates  the  difference  between  upstream 

13  * 

and  downstream  temperatures;  thus,  it  is  this  combination  which  is  used  to 
record  the  rate  of  change  of  upstream  and  downstream  temperatures  with  time, 
resulting  from  the  step  change  in  fluid  temperature  upstream.   This  differ- 
ential is  recorded  on  a  Minneapolis -Honeywell  "Brown"  Recorder  and  serves 
as  the  primary  data  for  heat  transfer  evaluation  of  the  matrix.   The  re- 
corder has  variable  chart  speeds  so  that  acceptable  cooling  curves  for 
maximum  slope  determination  can  be  generated.   This  instrument  has  a  span 
adjustment  which  permits  the  span  to  be  varied  continuosuly  from  the  narrow- 
est to  the  widest  span  desired  from  0  to  55  millivolts.   For  the  work  perform- 
ed herein,  the  recorder  was  precalibrated  for  a  0  -  3  millivolt  scale.   This 
unit  also  has  adjustable  supression  and  damping  adjustment  features.   The 
damping  adjustment  provides  a  filter  network  for  removing  a.c.  strays,  and 
provides  correct  damping  so  that  neither  undershoot  nor  overshoot  in  the 
response  curves  is  experienced. 

A  copper  constantan  thermocouple  was  inserted  into  the  duct  just  a- 

head  of  the  orifice  to  measure  T  ,  orifice  temperature.   This  temperature 

o 

was  read  in  millivolts  on  a  Rubicon  Company  Portable  Precision  Potentio- 
meter, converted  to  degrees  Fahrenheit  and  recorded  for  each  run. 

A  copper  constantan  thermocouple  grid  was  inserted  in  place  of  T 
to  determine  the  temperature  distribution  across  the  cross  section  during 
initial  check  out  of  the  apparatus.   This  thermocouple  arrangement  was 
part  of  a  special  holder  lined  with  balsa  wood  to  provide  a  smooth, 
continuous  flow  passage,  which  also  contained  an  impact  tube  used  for 
checking  the  velocity  profiles   (See  Figure  33  ). 

PRESSURE  MEASURING  SYSTEM:   (See  Figure  1) 

Pressure  taps  are  located  upstream  and  downstream  of  the  matrix  and  on 
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either  side  of  the  orifice.  Each  tap  is  connected  to  an  appropriate  mano- 
meter or  draft  gage  via  Imperial  Company  "poly-flo"  tubing.  Quick  closing 
valves  were  installed  at  various  positions  in  the  line  to  permit  isolation 
of  sections  and  as  a  safety  feature. 

The  following  instruments  were  utilized: 

(1)  Meriam  Instrument  Company,  0  -  50"  manometer 

(2)  Meriam  Instrument  Company,  -8"  to  +8"  manometer 

(3)  Ellison  Differential  Direct  Draft  Gages,  0  -  6",  0  -  8". 

(4)  Ellison  Inclined  Draft  Gages,  0  -  0.5",  0  -  6". 

(5)  E.  Vernon  Hill  and  Company  Type  "C"  Micromanometer, 
0  -  1.25". 

Any  one  of  these  instruments  or  a  combination  can  be  used  to  measure  the 
differential  pressures  of  the  orifice  or  matrix  or  the  required  static 
pressures  -  upstream  of  matrix  and  upstream  of  the  orifice.   Cross  checking 
of  the  various  instruments  assured  reliable  operation. 

FLUID  HEATER  SYSTEM:  (See  Figure  1) 

The  heater  section  is  comprised  of  28  nichrome  wire  heaters  of  0.0031" 
diameter.   The  heater  system  was  designed  to  elevate  the  air  temperature 
20  F  above  ambient  for  a  mass  rate  of  flow  of  1000  lb/hour.    The  nichrome 
material  was  selected  because  it  has  high  resistivity,  low  thermal  conduc- 
tivity and  specific  heat,  and  has  a  very  small  time  constant  so  that  it 
permits  one  to  approach  the  idealized  step  function.  [_l\ 

The  wire  heaters  are  wound  two  to  a  bakelite  frame,  1/32"  between  each 
wire,  with  50  and  52  wires  respectively  to  each  heater  on  a  frame.   These 
heaters  are  connected  in  parallel,  then  via  a  switch  to  a  variable  voltage 
bus.   The  heater  frame  and  its  electrical  connections  are  apparent  in 
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Figures  34  and  35. 

A  schematic  wiring  diagram  for  one  of  the  two  heater  complexes  is  shown 
below  |_1  J  .   A  total  of  14  frames  are  wired  in  this  fashion,  thereby  put- 
ting all  28  heaters  in  parallel,  permitting  two  to  be  switched  on  or  off 
at  a  time. 
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The  number  of  heaters  in  use  decreases  as  flow  rate  decreases.   The 
necessary  voltage  variations  are  obtainable  by  a  General  Radio  Company 
Type  W20HM  "Variac"  Autotransformer,  0  -  280  volt,  8  amp  load  from  a  240 
volt  50  -  60  cycle  line. 

The  recommended  procedure  is  to  decrease  the  number  of  heaters  insofar 
as  practical  while  maintaining  a  high  voltage.   When  there  are  only  four 
heaters  in  use,  voltage  reduction  becomes  necessary  to  obtain  the  required 
20   temperature  rise,  keeping  the  number  of  heaters  constant.   Meeting 
the  current  requirement  is  the  dictating  factor.   It  is  noteworthy  that  the 
heaters  should  never  be  energized  unless  sufficient  air  is  flowing,  other- 
wise the  wires  will  sag  and  short  circuiting  will  result. 
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MATRIX  HOLDER  AND  TEST  SECTION;   (See  Figures  31,  32,  and  35) 

The  matrix  holder  and  test  section  are  made  of  polyethylene  plastic. 
The  test  section  casing  holds  the  T  and  T,  thermocouple  grids  permanently 
mounted,  as  well  as  the  upstream  and  downstream  static  pressure  taps  required 
for  determining  friction  factor.  All  parts  were  machined  to  close  tolerance 
to  assure  a  snug  fit  for  the  matrix  holder  and  to  guarantee  good  alignment 
of  the  flow  channel  through  the  heater  and  test  sections.  The  flow  channel 
is  3  -  3/16"  x  3  -  3/16"  and  matrices  of  flow  lengths  up  to  3"  may  be 
tested.   As  mentioned  previously  under  Temperature  Measuring  System,  the  T 
thermocouple  grid  is  located  in  the  matrix  holder. 

Correct  positioning  of  the  test  matrices  in  the  flow  channel  is  en- 
sured by  using  styrofoam  plastic  inserts.   In  addition  to  the  aforemention- 
ed systems,  one  item  remains  to  be  considered  -  the  Inlet  Cone  and  Screen 
Straightener  Section.   The  inlet  cone  was  designed  using  Smith  and  Wang 
[l8j  as  a  guide.  Although  this  reference  pertained  to  contracting  cones  for 
circular  duct,  the  resultant  inlet  cone  for  this  rig  proved  very  satisfactory. 
To  assist  in  guaranteeing  a  uniform  velocity  at  the  throat,  a  stainless 
steel  60  mesh  screen  was  installed  (See  Figure  1).   Space  is  available  in 
the  screen  straightening  section  to  add  additional  screens  if  they  should 
be  required  for  flow  rates  beyond  the  range  of  this  investigation. 
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APPENDIX  C 
PERFORATED  NICKEL  GEOMETRIC  PROPERTIES 
The  material  used  in  the  perforated  plate  matrices  was  pure  nickel, 
electro-deposited  sheet  of  integral  structure  manufactured  by  Perforated 
Products,  Incorporated.   Types  160/40TV,  160/40Q,  50G,  125M  and  125P 
were  utilized  in  the  matrices  -  the  first  two  of  these  having  slotted 
openings  and  the  remaining  three  having  round  openings.   A  brief  descrip- 
tion of  each  of  the  aforementioned  types  as  specified  by  the  manufacturer 
is  set  forth  in  Table  C-l. 

Inasmuch  as  the  perforations  in  all  cases  were  conical,  the  use  of 
the  average  slot  length  or  average  hole  diameter  was  considered  a  propos. 
Close  examination  of  each  of  the  different  types  revealed  that  this 
average  value  was  somewhat  less  than  the  dimensions  specified  by  the  manu- 
facturer in  each  case.   In  addition  to  the  requirements  for  determining 
porosity  of  the  plate,  a  correction  for  the  increased  longitudinal  conduc- 
tion length  due  to  the  perforations  as  well  as  a  correction  to  the  solid 
conduction  cross-sectional  area  was  performed.   The  following  idealized 
patterns  for  each  type  used,  best  illustrate  these  corrections. 
160/40TV 
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Cross  Hatched  Area,  A   =  (.0063)  (.0.25)  =  .0001575  in. 
Slotted  Area,  A  ,,  in  Ach  =  (.0008)  (.0155)  =  .0000124  in. 
Plate  Porosity,  the  fraction  of  open  or  slotted  area, 

A 


7/       o 


Si 


lch 


Solidity,  s,  is  1-y  =  1  -  si 
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s  =  1 


ch 

.0000124 
.0001575 


=  0.9212 


.0095" 


te   If 


-  <= 


h- 


_^> 


0125" 
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.0063" 


Li 


L.  ,  conduction  path  length,  =  0.0125  +  0.0063  =  0.0188" 


l 


.0063 
.0188 


■  0.3351 


The  solid  cross  sectional  area  for  heat  conduction,  A   varies 

i 
along  the  path  length.   The  area,  A^  ,  selected  was  the  cross  sectional 

i 
area  between  the  slots  normal  to  the  flow  direction,  i.e., 

A   =  (0.0095)(0.0022)  =  0.0000209  in.2   ■ 
i 

Just  as  L  =  2  L  where  L  is  the  total  matrix  flow  length  and  L, 
is  the  flow  length  between  two  successive  perforations;  so  does  L,  =  2_.  L, 
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and  A.=  "22  A,  .  In  evaluating  A  ,  a  conduction  area  reduction  ratio  is 

X  /    cross  sectional  area  of  perforations; 

employed  and  is  defined  as  ?     =  n   ♦.•«„.,i — 7^7, 

r  *  '    cross  sectional  area 

therefore,  A  ,  is  the  solid  cross  sectional  area  multiplied  by  1  -£ 

i.e., 

A,  =  A  (1  -  ^  ) 


-  i  finsft   f 1   (.0155) (.0022)1  „  Q   6102 
-1.6058   ^1  -  (.025)(.o022)  J 


in 


-  0.00424  ft, 


160/40Q 


A  ,_  -  (.0063)(.025)  =  .0001575  in, 

ch 


A  1 

si 

= 

(.0019)(.0161) 

-  .00003059 

in 

1/ 

■ 

.1942 

s 

ss 

.8058 

L 

33 

0.3351, 

same  as  160/40TV 

Lk 
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■ 

.644 
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■ 
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./>  - . 

002887  ft.2 
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125M 


k-  .0079" 


a  =  .00366" 


.0079  sin  60° 

=. 00634" 


A   =  (.00684) (.0079)  =  .00005404  in. 
ch 


lhole 


(.00366)2  =  .00001052  in. 


7/ 


We   =  .00001052   -  0.19466 


ch 
s   =  .8053 


.00005404 


I*.  .0079"  —J 


T 


Li  =  .0079  sin  60° 
=.00684" 


) 


97 


The  equivalent  conduction  length,  L.  ,  is  computed  in  the  case  of 

Ki 

round  perforations  assuming  the  path  to  be  as  shown  above,  i.e.,   2  TT  r^ 

3 

where  r   is  the  radial  distance  equal  to  one-half  the  center  to  center 

distance.   Thus 

L.   =  2TY   (.0079)=  .00827" 
ki         3      2 

and 

Li 

L   =  .00684  =  .82709 
i   .00827 

A  is  determined  as  in  the  case  of  the  slotted  perforations, 
k 

p       =  diam.  of  hole  x  plate  thickness 

center  to  center  distance  x  plate  thickness 

-  (.00366)00016)  =   .4633 
(.0079)00016) 

Therefore,     A  -  A  (1  -  $       )  =  1.1679in.2  (1  -  .4633) 
k    s 


=  .6268  in.2  =  .00435  ft.2 


125P 
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Ach  =  (•00684)(-0°79)  =  .00005404  in2 


50G 


'hole 


(.00291)   -  .00000665 


Z'  =  0.12307 
s  =  0.8769 


L.  .  =  .00827" 

Li  =  .82709 

Lki 


j>   =  0.3684 
A  =  A  (1  -  f     )  =  .00640  ft2 


|*_..0197" 


d  =  .0138" 


1^  =  .0197  sin  60° 
=  .01706" 
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Ach  -  (.01706) (.0197)  -  .00033608  in.2 

2  2 

Ahole  =  %—     (-0138)  ■  .00014957  in. 
4 

Z'  =  .44504 
s  -  .55496 


L.   =  .02063" 
ki 


li     =  .82697 
Ki 

$     -  .7005 
A.  -  Ag  (1  -  f      )  =  .00243  ft.2 
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APPENDIX  D 
SAMPLE  CALCULATIONS 
The  following  sample  calculations  are  based  on  data  obtained  from 
Run  #1  for  the  modified  parallel  plate  matrix  constructed  from  the  160/40Q 
perforated  nickel  material.   The  basic  geometric  parameters  are  illus- 
trated in  Appendix  A. 

Prior  to  the  reduction  of  data,  several  fixed  parameters  must  be 
evaluated. 

DIMENSIONS  OF  FINNED  SHEETS:   3.9375"  x  2.00"  x  .0016" 
(BEFORE  FORMING) 

DIMENSIONS  OF  SPLITTER  PLATES:   3.1875"  x  2.00"  x  0016" 

NUMBER  OF  FINNED  SHEETS:   102 

NUMBER  OF  SPLITTER  PLATES:   103 

WEIGHT  OF  MATRIX:   311.0979gms/453.6gm/lb  -  0.68584  lb. 

MATERIAL  CONSTANTS: 

k  =38.7  Btu/hr  -   F  -  ft 
s 

c  =  0.1065  Btu/lb  -   F 
s 


?  =  0.321  lb/in3  =  554.7  lb/ft3 


PLANE  SURFACE  AREA,  A*,    -   (102) (2) (2) (3. 9375  +  3.1875) 

144 

20.1875  ft2 

HEAT  TRANSFER  AREA,  A,   =  PLANE  SURFACE  AREA  x  SOLIDITY 

A  =  A*s 

=  (20. 1875) (.8058) 


A  =  16.2667  ft2 


FRONTAL  AREA,  A  ,  =  3.141"  x  3.1875"  =  10. ,012  in2 

fr 


Afr  =  0.06953  ft 
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It  has  been  shown  in  Appendix  A;  however,  that  this  can  be  alter- 
nately expressed  as: 

D  =  4_e_ 

Inasmuch  as  the  perforations  have  no  effect  on  the  hydraulic  diameter 
the  value  of  Qi  rather  than  Q     is  used. 


DR  -  (4)   (0.8834)  = 
(1742.1) 


0.002028  ft 


Now,  we  are  in  a  position  to  reduce  the  data.   The  following  presenta- 
tion outlines  the  steps  required;  however,  as  far  as  practicable  the 
actual  evaluation  is  performed  by  the  digital  computer  program  set  forth 
in  Appendix  E. 


RECORDED  DATA: 


d  =  2.310" 
o 

AP  =  6.69  in  H20 

A  P  -  10.15  in  HO 

m  2 

P  ■  10.65  in  HO 
s  2 


Po  -  20.80  in  HO 


T  =  69.0  °F 
o 

Patm-  30-140  inHg' 
Chart  Speed  =  4  in/sec. 

(3  -  d  /d  =  2.310  =  0.75 
0     3.08 


DETERMINATION  OF  MASS  RATE  OF  FLOW:       (See  Appendix  A) 

2 


m   -     589.81  Fa  Y  d, 


fP? 


AP      x    P. 
o        _12 

T 
o 


where:      C  =  C, 


+  AC   f    10 

L     NR(d   ) 


from  reference 


M 


P     -   (P  ,      -  P  /13.6)   0.4912  I  T     -  t     (°F)  +  459.7 
3  atm         o  o  o 

Y     ^1.0     Fig.   40b  of  Reference    [l6] 
F       =   1.0   Fig.    30  of  Reference    [l6] 
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For  g  =0.75  and  d  »  3.08",   C  -  0.60691  and  AC  =  0.03839,  respective- 


[■ 


P,  =  30.140  -  20.80 
13.6 


x   0.4912    -  14.0537  psi 


SOLID  MATRIX  CROSS  SECTIONAL  AREA,  Ag, 

As  =  102  x  3.9375"  x  .0016"  +  103  x  3.1875"  x  .0016" 


=  1.1679  in2 


Ag  =  0.00811  ft' 


FREE  FLOW  AREA,  A.,  -  FRONTAL  AREA  -  SOLID  CROSS  SECTIONAL  AREA 


A„  =  A,  -  AQ 
c    fr    s 


Ac  =  0.06142  ft' 


CONDUCTION  AREA  CORRECTED  FOR  EFFECTS  OF  PERFORATIONS,  A 

K 

A^  ■  A  (1  -  p   )  where  *     is  the  conduction  area  reduction 
8  ratio  (See  Appendix  C). 

-  0.00811  (1-0.644) 


Ak  =  0.00288  ft' 


POROSITY,  p,  =  MATRIX  FLOW  VOID  VOLUME  -  Vm  -  Vs  =  1  -  Vs 

MATRIX  VOLUME  Vm  Vm 


0.8834 


p  =  1  -  0.00135 
0.01159 


COMPACTNESS,  (3   =   TOTAL  HEAT  TRANSFER  AREA 

MATRIX  VOLUME 


A 
Vm 


16.2667  = 


1403.75  ft2/ft3 


0.01159 

Note  that  (3   includes  the  effect  of  area  reduction  due  to  perfora- 
tions. 

(3  for  an  unperf orated  surface  would  therefore  simply  be  the  plane 
surface  area  divided  by  the  matrix  volume. 
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3  -a* 

s    Vm 


20.1875  =  1742.1  ft2/ft3 
0.01159 


THE  HYDRAULIC  DIAMETER,  D  ,  =  4rL  = 

"      H'  h 


4   FLOW  CROSS   SECTIONAL  AREA 
WETTED   PERIMETER 


T     =   69.0  +  459.7  =  528.7     °R 
o 

m    -     (589.81   (2.310)2 

!  1  -  (0.75)4 
m«  1605.242  C 
assume  C  =  0.6000 

then  rr\   =  963.145  lb/hr 
C  -  CQ  +  L  C  J   lO4, 

L  Nr<"o). 


\  6.69  x  14. C 
\       528. 


0537 

7 


where   N, 


R(d0)  =  (4)  (12)  m 

(3.1416)(3.08)/a 

=  4.961  20. 


u  is  evaluated  at  temperature  t  by  means  of  the  linear  approximation 
IX   =  0.0395  +  0.64167  x  10'4  tQ  =  0.0395  +  0.64167  x  10"4  (69.0) 

=  0.04392  lb/ft  -  hr 
therefore, 

C  =  0.60691  +  0.03839 
C  =  0.6104 
This  value  of  C  is  compared  with  the  assumed  value  of  C  and  successive 
iterations  are  performed  until  the  two  values  are  equal.   The  iterative 
process  shows  that  a  value  of  C  =  0.6172  satisfies  this  requirement  and 
thus 


(10^) (0.04392) 
(4. 961) (963. 145) 


m  =   990.756  lb/hr 
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MATRIX  REYNOLDS  NUMBERS: 


Two  values  for  Reynolds  number  are  required,  one  for  the  isothermal 
flow  friction  factor  evaluation  and  second  for  the  heat  transfer  evalua- 
tion. The  two  values  differ  only  in  so  far  as  the  difference  in  the 
absolute  viscosities.  JJin     is  evaluated  at  orifice  temperature  (t  ) 
whereas  JUH  (heat  transfer)  is  evaluated  at  average  bulk  fluid  tempera- 


ture which  is  assumed  to  be  at  t  +  10  F. 

o        o 


NRf  =  £LjL 


NRH  =  ^^ 


yH 


(4)   12   (990.756)  = 
(16.2667) (.04392) 


924.54 


=  NRf^£ 
/Jh 


=  924.54  (.04392)  = 
(.04456) 


911,23  | 


CONDUCTION  PARAMETER: 


therefore 


\  =   ks  Ak   -  (38.7)(0»00288) 

LAc 
P 


=  0.00282 


(2/12) (990. 756) (.24) 


k 


Xjt- 


L     =  0.3351  from  Appendix  C. 


\v. 


(.00282) (.3351)  = 


0.000944 


MAXIMUM  SLOPE: 

The  slope  of  the  cooling  curve  was  determined  by  the  method  pres- 
cribed in  Appendix  A  from  Figure  36.   Slopes  were  determined  for  each  of 
the  curves  and  the  average  value  was  used. 


MAXIMUM  SLOPE  =  C_ 


[*] 


max 


x  SLOPE  x  CHART 
SPEED 


C  -  W  c 
s    s  s 


(0.68584  lb) (0.1065  Btu/lb  °F)  -  0.07304 


C  =  mc   =   (990.756  lb/hr)(0.24  Btu/lb  *F)  __1_ 


hr   =  0.06605 
3600    sec 
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0.63507 


SLOPE  =  0.14357   (See  Figure  36) 

whence, 

MAXIMUM  SLOPE  =  (.07304)   (0.14357)  (4)  = 
(.06605) 

COLBURN  j  FACTOR: 

Enter  Table  XI  with  the  values  of  V\   and  MAXIMUM  SLOPE  to  deter- 

k 

mine  N_   ■  4.05.   Linear  interpolation  from  this  table  is  sufficient. 

j  ■  N.ju   Ac     NpR  2/3   from  Appendix  A, 

"K 

2/3 

N        -  0.8017  -  0.82353  x  10    t0  using  a  linear  approximation; 
PR 

therefore 


0.01217 


j  =  (4.05)(0.06142)   (.7960)  = 
(16.2667) 

FANNING  FRICTION  FACTOR: 

From  Appendix  A,  the  friction  factor  is: 

Patm  =  30'140  in  Hg  x  0.4912  £si    =»  14.8048  psi 

in  Hg 

P  -  10.65  in  HO  x  5.204  psf     =  55.42  psf  =  0.3849  psi 
S  Z         in  H20 

P  -  10.15  in  H90  x  5.204  psf      =  52.82  psf  =  0.3668  psi 

m  *  in  H20 

P,  =  P  fc  -  Pc  =  14.8048  -  0.3849  =  14.4199  psi 
1    atm    s  r 

P2  =  PL-  APm  =  14.4199  -  0.3668  =  14.0531  psi 


p 

mean 

=  Px     +     P             -  14.2365  psi 
2 

p 

mean 

■   ttmcan*    (144)            =     (14.2365) (144) 
R(t     +  459.7)                <53.35)(528.7) 
o 

G  =    m    =  990.756   =  16130. 841b/hr  -  ft2  =  4.481  lb/sec  -  ft2 

A 0.06142 

c 
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r 

H    ■  0.000507  =  0.00304 
7"     2/12 


K  =  0.48  and  K   =  -0.33  from  Figure  38. 
c  e 

(  1  +  p  )  =»  1  +  (0.8834)2  -  1.7804 


f  = 


(2)  (32.2)(0.07268)(52.82)  -   (0.48  -  0.33)  -  .3668 
(4.481)z  14.2365 

(1.7804)     (0.00304) 

(12.3126  -  0.150  -  0.04587) (0.00304) 


f  =  0.03686 


The  thermal  properties,  uand  NpR  '   ,  were  obtained  by  linear 
interpolation  of  data  from  reference  ^6  J  .  The  values  of  k  and  c 
were  obtained  directly  from  references  [5  1  and  [4  I  for  nickel  and 
stainless  steel  respectively. 
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APPENDIX  E 
FORTRAN  60  DATA  REDUCTION  COMPUTER  PROGRAM  FOR  A  CDC  1604  DIGITAL  COMPUTER 

A  computer  program  was  formulated  in  order  that  as  much  data  reduction 
as  practicable  would  be  performed  by  the  digital  computer.   (See  Table  E- 
1). 

It  was  hoped  that  this  program  together  with  a  subroutine  would  produce 
Colburn  j  factor  directly.   An  unsuccessful  attempt  was  made  to  fill  out 
table  values  of  N™   (TABLE  XI)  by  curve  fitting  using  Aitken's  iterative 
method  for  fitting  an  nth  degree  Lagrange  polynomial  to  the  n  +  1  given 
points  -  allowing  for  unequal  spacing.   As  a  result  of  reference  [8}  ,  this 
approach  was  abandoned;  however,  the  computer  program  presented  is  still 
helpful  in  that  values  of  conduction  parameter  and  maximum  slope  are  end 
results. 

Entering  TABLE  XI  with  these  parameters  yields  a  value  of  N  ,  which 

when  multiplied  by  CI  (CI  ■  Np  '   »  A  /A)  produces  the  Colburn  j  factor. 

2/3 
A  linear  approximation  for  values  of  viscosity  and  N       as  functions 

of  temperature  is  included.   An  added  convenience  in  the  computer  program 

is  the  evaluation  of  the  orifice  constant,  K.   Using  \_15 J  ,  an  assumed 

value  for  C  (a  slowly  varying  function)  was  made.  The  correct  value  for 

C  was  then  evaluated  by  an  iterative  process  by  the  computer.   With  the 

value  of  C  determined,  K  (a  rapidly  varying  function)  was  found  from  the 


following  simple  relationship:    r^  -   C 


where   fi   is  the  ratio  of  the  orifice  diameter  to  the  inside  pipe  dia- 
meter. 

With  the  computer  program  presented,  an  unlimited  number  of  matrices 
can  be  reduced  by  simply  changing  the  first  card  of  the  program.    (L  = 
number  of  matrices  for  which  data  is  to  be  reduced). 
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The  following  is  a  list  of  the  input  data  for  the  FORTRAN  60  program: 
(See  sample  input  data  in  TABLE  E  -  2) 
CARD  NO. 

30        XML       Matrix  length,  ft 

o 
AS        Solid  conduction  area,  ft 

RH        Hydraulic  radius,  ft 

POR       Porosity 

2 
AFL       Free  flow  area,  ft 

2 
AHT       Matrix  heat  transfer  area,  ft 

WM  Matrix  weight,  lbs. 

40        EK  Exit  flow  coefficient 

CK  Entrance  flow  coefficient 

SK  Solid  thermal  conductivity,  Btu/hr  ft  °  F 

CM  Matrix  material  specific  heat,  Btu/lb  °F 

XLR       Ratio  of  L/I*,   (Flow  length/Equivalent  conduc- 
tion length) 

N  Number  of  runs 

10        DO  Diameter  of  orifice  plate  opening,  inches 

DELPO  Pressure  drop  across  orifice,  inches  H  0 

PO  Static  pressure  upstream  of  orifice,  inches  H20 

ATMP  Atmospheric  pressure,  inches  Hg 

BETA  Ratio  of  orifice  diameter  to  inside  pipe  diameter 

10+2    CS  Chart  speed,  sec/inch 

Pressure  drop  across  matrix,  inches  H_0 


DELII 


HS        Static  pressure  at  matrix  inlet,  inches  H  0 
SLO       Slope  of  generalized  cooling  curve,  (inches)"* 
TEMPO     Fluid  temperature  at  orifice  (t  ),   °F 
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The  output  from 
EMDOT 
RNUMP 
RNUMH 
CAPFL 
CAPS 
CONDPAR 
SLOMAX 
CONPARE 
G 

RNUMFR 
FFR 

ORIFICE  K 
PRANDTL  NO. 
CI 
UPR 

ESTD 
STDH 


the  program  is  as  follows: 

Mass  rate  of  flow  (m),  lb/hr 

Reynolds  number  for  the  pipe 

Reynolds  number  for  heat  transfer 

Thermal  capacitance  of  the  fluid  (C),  Btu/sec   F 

Thermal  capacitance  of  the  solid  (C~_),  Btu/  °F 

Conduction  parameter  (  A.  ) 

Maximum  slope  of  generalized  cooling  curve 

Equivalent  conduction  parameter  (  A»K    ) 

Mass  flow  velocity  (m/A  ),  lb/sec  -  ft 

Reynolds  number  for  friction 

Fanning  friction  factor 


Orifice  constant,  K 

>er)' 
2/3 


2/3   (Prandtl  number) 2/3  N„  2/3 

>   Pr 


Parameter  =  N  "'  "  x  A„/A 
Pr       c 

Fluid  viscosity,  yW   ,  lb/hr-ft 

(evaluated  at  fluid  temperature  at  orifice) 

Flow  friction  power  per  unit  area  (E   ),  HP/ft2 

Parameter  for  evaluating  heat  transfer  power  per  unit 
area  (heTn),  Btu/hr  ft2  °F.   (STDH  x  NTi  =  hqTn) 


lSTDJ 


'Tu   "STDJ 
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PROGRAM  SING3LO 
L=l 
K=l 
30  READ  300»XML»AS»RH»POR»AFL»AHT>WM 

300  FORMAT (  7F  10. 5) 

40  READ  400,EK,CK,SK,C-1,XLR,N 
400  FORMAT ( 5F10.5. 12) 
M=l 

10  READ  10G,DO,DELPG,PO,ATMP,3ETA 

100  FORMAT(5F10.5) 

READ  101,CS»DELH  ,HS * SLO , TEMPO 

101  FORMAT ( 5F 10. 5) 

P3=( ATMP-(PO/13.6) )*.4912 
T0=TEMPO+459.7 
UFR =0.0 395+0. 64167 E-4*TEMP0 
T3=TEMPO+10.0 
UH=0.0 395+0. 64167E-4*TB 
WRITE  OUTPUT  TAPE  4,9 
■  9  FORMAT ( 1H1) 

WRITE  OUTPUT  TAPE  4,4,M 
4  FORMAT (9H  RUN  NO.  12///) 

11  IF(ABSF(BETA-.75)-.000  001)27»27»12 

12  IF(ABSF(BETA-.64)-.000001)26»26»13 

13  I F  (  ABSF  (  BET  A-.  50  j  -  .000  001 )  25'»  2  5  » 14 

14  IF(lABSF(BETA-.40)-.000001)24.24,15 

15  I F< ABSF ( BET A-. 25) -.000001) 23 » 23 » 16 

16  IF(ABS  F(l3ETA-. 15)-. OOOOol  122,22,17 
17.  GO  TO  8  0 

21  X2=0. 59225 
DELX=0. 01916 
ASSUM=0.60  0  0 
GO  TO  2.8 

22  X2=0. 59171 
DELX=0. 01691 
ASSUM=0.6000 
GO  TO  28 

23  X2=0. 59184 
DELX=0. 01622 
ASSUM=0.6000 

GO  TO  28  ■' 

24  X2=0. 59448 
DELX=0. 01867 
ASSUM=0.6000 

■  GO    TO    28 


TABLE  E-l 
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2  5    X2 =0*598  50 
DELX=0. 02125 

ASSUM=0.6000    '  t 
GO    TO    28 

26  X2=C. 60575 
DELX=0. 02595 
AC,SUM=0.60C0 
GO    TO    28 

27  X2=0. 60691 
DELX=0. 03839 
ASSUM= 0.600  0 

23    COR=ASSUM 
NC  =  0 
B1=1../SQRTF(  l.-BETA**4) 

31  NC=NC+1 

IF (NC-4) 29*29.60 
2  9    EMDOT=5  89.81*DO*DO*SORTF( DELPO*P3/ TO ) *COR*dl 
RNUMO= ( 4.961*EMDOT/UFR ) /3ETA 
COR A=X2+DELX -SORT F( 10000. /RNUMO ) 
WRITE    OUTPUT    TAPE    4, 8, CORA 
8    FORMAT(6H    COR A= F 10 . 5 / / ) 

32  COR=CORA 
GO    TO    31 

60    PR  =  0.8017-.82353E-4*TEivPO 
C1=PR*AFL/AHT 
CP=0.24 

RNUMP=5.09*EMDOT/UFR 
RNUMH={4.0*XML*EMDOT )/ (AHT*UH) 
RNUMFR=RNUMH*UH/UFR 
CAPFL=EMDOT*CP/3600.0 
CAPS=WM*CM 

SLOMAX= (CAPS*SLO) / (CAPFL*CS ) 
CONDPAR=(SK*AS) / < XML*EMDOT*CP ) 
COIMPARE=CONDPAR*XLR 
B=CK+EK 
C=1.0+POR**2 
G=EMDOT/{AFL*3600.0) 
DELP=DELH*. 03613 
PS=HS*. 03613      • 
PA=ATMP*.4912 
P1=PA-PS 
P2=P1-DELP 
PM=(Pl+P2)/2. 
RHOM=( PM*144.0 ) / ( 53. 3* ( TEMPO  +  460 . 0  )  ) 


TABLE  E-l 
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FF1={ (64.4*RHOM*DELP*144.0 )/ (G*G) ) 
FF2=ODELP/PM 
FF3=FF1-FF2-B 
FFR=FF3*RH/XML 
C0NST=C0R*B1 

STDH=(0.02195*RNUMH*C1  )/  (4.*RH) 

ESTD=1.  ll£-7*<  1.0/  (4.*RH)  ) **3*FFR* ( RNUMFR/1000.  )***3 
WRITE    OUTPUT    TAPE    4 , 7 , STDH * ESTD 
7    F0RMAT(19H    CONSTANT    FOR    hSTD=F 10 . 5 , 3X , 6H    ESTD=F10 . 5/// ) 

WRITE    OUTPUT    TAPE    4> 1 »EMDOT »RNUMP>RNUMH*CAPFL»CAPS 
1    FORMAT (7H    EMD0T=F10 . 5 , 3X » 7H    RNUMP=E20.8 »3X»7H    RNUMH=F10. 5 » 3X • 
17H    CAPFL=F1C.5»3X»6H    CAPS  =  F10. 5/// ) 
WRITE    OUTPUT     TA^E    4* 2 » CONDPAR , SLOMAX *  COMPARE 
2-F0RMAT(9H    C0NDPAR=F10 . 5 » 3X »8H    SLOMAX  =  F10.5»3X >9H    CONPARE=F 10  .  5/ //  ) 

WRITE    OUTPUT    TAPE    4 , 3 ,G » RNUMFR , FFR , CONST 
30FORMAT(3H    G=F10 . 5 » 3X » 8H    RNUMFR  =  F1Q, 5>3X » 5H    FFR= F10  .  5  ,  3X  ,  1  lri    ORIFIC 
IE    K=Fl0.5///) 
WRITE    OUTPUT     TAPE    4 , 5  ,  PR » C 1 , UFR 
5    FORMAT (16H    PRANDTL    NO. 2/3=F10. 5 » 3X »4H    C l=r 10 . 5  » ^X , in    UFR=F10.5/// ) 
M  =  M+1 

IF(N-M) 50,10,10 
50    K  =  !<+1 

IF(L-I<).20,30,30 
8  0    PRINT    81 
81    FORMAT (30H    ORIFICE    CONSTANT    ERROR    EXISTS) 

GO    TO    20 
20    STOP 
END 
END  < 
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SAMPLE  INPUT  DATA  FOR  PROGRAM  SINGBLO 


0,1667 

.00640 

.00049 

.8542 

.05939     17.70242 

-.32 

.51 

38.7 

.1065 

.82709    11 

2.310 

6.78 

21.5 

30.140 

.75 

.25 

9.30 

12.2 

.17768 

68.5 

2.310 

4.40 

15.45 

30.120 

.75 

.25 

7.15 

8.30 

.15259 

71.0 

2.310 

2.74 

10.95 

30.120 

.75 

.25 

5.15 

5.80 

.12882 

73.0 

1.971 

3.82 

7.10 

30.120 

•  64 

.50 

3.40 

3.50 

.21040 

72.0 

1.540 

6.89 

4.50 

30.120 

.50 

.50 

2.33 

2.06 

.17379 

72.7 

1.540 

2.91 

2.60 

30.110 

.50 

.50 

1.42 

0.95   • 

.14039 

72.7 

1.232 

2.75  . 

1.40 

30.110 

'.40 

.50 

0.80 

0.43 

.09232 

72.5 

0.775 

6.00  • 

0.68 

30.110 

.25 

1.00 

0.41 

0.18 

.10170 

69.0 

0.775 

5.98 

0.67 

30.110 

.25 

1.00 

0.40 

0.178 

.10316  - 

.  69.0 

0.775 

1.78 

0.35 

30.110 

.25 

2.00 

0*22 

0.075 

.10321 

70.0 

0.462 

2.60 

0.15 

30.110 

.15 

4.00 

0*080 

0.025 

.07513 

70.0 

•  4844 
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